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ABSTRACT 

We analyze five years of PLANET photometry of microlensing events toward the Galactic 
bulge to search for the short-duration deviations from single lens light curves that are indicative 
of the presence of planetary companions to the primary microlenses. Using strict event selection 
criteria, we construct a well defined sample of 43 intensively monitored events. We search 
for planetary perturbations in these events over a densely sampled region of parameter space 
spanning two decades in mass ratio and projected separation, but find no viable planetary 
candidates. By combining the detection efRciencies of the events, we find that, at 95% 
confidence, less than 25% of our primary lenses have companions with mass ratio q = 10~^ and 
separations in the lensing zone, [0.6 — 1.6]6'e, where 6^ is the Einstein ring radius. Using a 
model of the mass, velocity and spatial distribution of bulge lenses, we infer that the majority 
of our lenses are likely M dwarfs in the Galactic bulge. We conclude that < 33% of M-dwarfs in 
the Galactic bulge have companions with mass nip = Mjup between 1.5 and 4 AU, and < 45% 
have companions with nip — 3Mjup between 1 and 7 AU, the first significant limits on planetary 
companions to M-dwarfs. We consider the effects of the finite size of the source stars and 
changing our detection criterion, but find that these do not alter our conclusions substantially. 

Subject headings: gravitational lensing, planetary systems 
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Introduction 



The discovery in 1995 of a massive planet orbiting 51 Peg (Mayor & Queloz 1995), followed by the 



discovery of many more planets orbiting nearby dwarf stars using the same radial velocity technique (Marcy, 



Cochran fc Mayor 2000 and references therein) has focussed both public and scientific attention on the 
search for extrasolar planets and the experimental and theoretical progress being made in developing other 
viable detection techniques. 

Due to their small mass and size, extrasolar planets are difficult to find. Proposed detection methods 
can be subdivided into direct and indirect techniques. Direct methods rely on the detection of the reflected 
light of the parent star, and are exceedingly challenging due to the extremely small flux expected from 



the planet, which is overwhelmed by stray light from the star itself ( Angel fc Woolf 1997 ). Some direct 



imaging searches have already been performed (Boden et al. 1998), but the future of this method lies in the 



construction and launching of space-based instrumentation (Woolf & Angel 1998) 



Astrometric, radial velocity, and occultation measurements can be used to detect the presence of a 
planet indirectly. Astrometric detection relies on the measurement of the positional wobble of the stellar 
centroid caused by the motion of the star around the center of mass of the planet-star system and yields the 
mass ratio and orbital parameters of the planet-star system. Many attempts to find extrasolar planets in this 
way have been made, but the measurements are difficult and the detections remain controversial; planned 
space-based missions astrometric missions such as the Full-Sky Astrometric Mapping Explorer (FAME), the 
Space Interferometry Mission (SIM), and the Global Astrometric Interferometer for Astrophysics (GAIA) 
are expected to be substantially more successful. Occultation methods use very accurate photometry of 
the parent star to detect the small decrease in flux (^ 1%) caused by a planet transiting the face of 



the star (Borucki & Summers 1984; Hale & Doyle 1994). Many occultation searches are currently being 
conducted (Deeg et al. 1998; Brown fc Charbonneau 200(]| ), with important new limits being placed on 
planetary companions in 47 Tuc ( iGiUiland et al. 200o| )~ Recently, one of the extrasolar planets detected 
via radial velocity surveys was also found to transit its parent star, yielding a measurement of the mass, 
radius, and density of the companion (Charbonneau et al. 2000| ; [Henry et al. 2000|). Spaced-ba sed missions 



are being planned to increase the sensitivity to low- mass planets (COROT, Deleuil et al. 1997 ; KEPLER 



Borucki et al. 1997). By far the most successful indirect method for discovering planets has been the 



Doppler technique, which employs precise radial velocity measurements of nearby stars to detect Doppler 
shifts caused by orbiting planets. Several teams have monitored nearby stars with the aim of detecting the 
Doppler signal of orbiting planets ( McMillan et al. 1993| ; Mayor fc Queloz 1995; Butler et al. 1996; Cochran 
et al. 1997 ; Noyes et al. 1997 ; Vogt et al. 200C ). To date these groups combined have discovered over 50 
extrasolar planets, with new planetary companions being announced every few months. Several exciting 
discoveries using the radial velocity technique include the flrst detection of extrasolar planetary systems 



( [Butler et al. 1999| ; [Marcy et al. 2001a ; Marcy et al. 20011 ; Fischer et al. 2002[ ) and the detection of planets 
with masses below that of Saturn (Marcy, Butler, & Vogt 2000). 



These detection techniques are complementary to one another both in terms of their sensitivity to 
planetary mass and orbital separations and the speciflc physical quantities of the planetary system that they 
measure. All share two distinct advantages: the experiments are repeatable and, due to their reliance on 
flux measurements of the parent star or the planet itself, they are sensitive to stars in the solar neighborhood 
where follow-up studies can be easily pursued. For example, spectroscopic follow-up studies may enable 
the detection of molecules commonly thought to be indicative of life, such as water, carbon dioxide, and 



ozone ( Woolf fc Angel 1998 ). This advantage is linked to a common drawback: most of the searches can be 
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conducted only on a limited number of nearby stars, and are thus unable to address questions about the 
nature of planetary systems beyond the immediate solar neighborhood. In addition, most of the methods 
(astrometry, radial velocity and occultation) can only probe companions with orbital periods smaller than 
the duration of the experiment. Furthermore, most are fundamentally restricted to massive planets, for 
example, radial velocity searches probably have an ultimate limit of 1 m due to random velocity 



variations intrinsic to the parent stars (Saar, Butler, fc Marcy 1998|). Of these methods, only space-based 



interferometric imaging and transit searches are expected to be sensitive to Earth-mass planets. 

Microlensing is a relatively new method of detecting extrasolar planets that overcomes many of 
these difficulties. Galactic microlensing occurs when a massive, compact object (the lens) passes near the 
observer's line of sight to a more distant star (the source). If the observer, lens, and source are perfectly 
aligned, then the lens images the source into a ring, called the Einstein ring, which has angular radius 



7E = 
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where M is the mass of the lens, is defined by, 
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and Dl and -Ds are the distances to the lens and source, respectively. The lens-source relative parallax is 
then TTrei = AU/Drei- Note that 9e corresponds to a physical distance at the lens of 



r. = 0.D.^2AV^^J . (3) 

If the lens is not perfectly aligned with the line of sight to the source, then the lens splits the source into 
two images. The separation of these images is ~ 26'e and hence unresolvable. However, the source is also 
magnified by the lens, by an amount that depends on the angular separation between the lens and source in 
units of 0E- Since the observer, lens, and source are all in relative motion, this magnification is a function 
of time: a 'microlensing event.' The characteristic time scale for such an event is 



.e^^=20da,s(^^J , (4) 

where fi^ci is the lens-lens relative proper motion, which we have assumed to be typical of events toward the 
Galactic bulge, ^rci = 25 kms~"^ kpc~^. 

If the primary lens has a planetary companion, and the position of this companion happens to be 
near the path of one of the two images created during the primary event, then the planet will perturb the 
light from this image, creating a deviation from the primary light curve. The duration tp of the deviation 
is roughly the time it takes the source to cross the Einstein ring of the planet, 9p. From equation (|^), 
9p = (r7ip/Af)^/^6'E, where nip is the mass of the planet. Therefore, from equation tp — {rUp/My/'^tE, 
or 

tp = V^tE (5) 

where q = rUp/M is the mass ratio of the system. For a Jupiter/Sun mass ratio (g ~ 10~^), the perturbation 
time scale is C'(day). Since the perturbation time scale is considerably less than ^e, the majority of the 
light curve will be indistinguishable from a single lens. Hence the signature of a planet orbiting the primary 
lens is a short-duration deviation imposed on an otherwise normal single lens curve. 
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Because microlensing relies on the mass (and not light) of the system, planets can be searched for 
around stars with distances of many kiloparsecs. Also, the sensitivity can, in principle, be extended down 



to Earth-mass planets (Bennett fc Rhie 1996). Finally, orbital separations of many AU can be probed 



immediately, without having to wait for a full orbital period. The primary disadvantages of microlensing 
searches for planets are that the measurements are not repeatable and there is little hope for follow-up 
study of discovered planetary systems. 
I I 

Mao & Paczyhski (1991) first suggested that microlensing might be used to find extrasolar planets. 

Their ideas were expanded upon by Gould fc Loeb (1992 ), who in particular noted that if all stars had 
Jupiter-mass planets at projected separations oi ^ rE, then ^ 20% of all microlensing events should exhibit 
planetary perturbations and that the detection probability will be highest for planets with projected 
separations lying within [0.6 — 1.6]0e of the primary, the "lensing zone." Since these two seminal papers, the 
theoretical basis of planetary microlensing has developed rapidly. Numerous authors have studied detection 



Bennett 



probabilities and observing strategies incorporating a variety of new effects (Bolatto & Falco 199S; 
fe Rhie 1996|; |Peale 1997|; ^ackett 1997|; [Griest fc Safizadeh 1998| ; |Gaudi, Naber , fc Sackett 1998|; [Pi Stefano 



fc Scalzo 1999a| ,|^; [Vermaak 2000| |Han fc Kim 200 1| ; [Peale 2001 1 ). Notably, [Bennett 
that the detection probability for Earth-mass planets could be appreciable (~ 2%) 



Rhie (1996) found 



and 



Griest & Safizadeh 



(1998) found that for high magnification events the detection probability can be nearly 100% for Jovian 



planets in the lensing zone, paudi fc Gould (1997 ), Gaudi (1998 ) and Gaudi fc Sackett (200C ) all discussed 
extracting information from observed microlensing events. In particular, Gaudi fc Sackett (2000| ) developed 
a method to calculate the detection efficiency of observed datasets to planetary companions; this method is 
employed extensively here. Planetary microlensing has been placed in the global context of binary lensing 
by Dominik (1999bD , and studied via perturbative analysis by Bozza (1999| , 20004 [b|)- 



On the observational front, progress has been somewhat slower. This is primarily because the survey 
collaborations that discover microlensing events toward the Galactic bulge, EROS ( Derue et al. 199"s| ), 
MACHO (lAlcock et al. 1997a|) , and OGLE ( [Udalski et al. 2000| ), have sampling periods that are of order 
or smaller than the planetary perturbation time scale, tp. However, soon after these searches commenced, 



these collaborations developed the capability to recognize microlensing events in real time (Alcock et al 



1996; Udalski et al. 1994), thus allowing publically available alerts of ongoing events. In response to this 



potential, several "follow-up" collaborations were formed: OMAN ( [Pratt et al. 1996| ; [Alcock et al. 1997b[ ), 
PLANET ([Albrow et al. 1998[) and MPS ([Rhie et al. 1999a[), with the express purpose of intensively 



monitoring alerted events to search for deviations from the standard point-source point-lens (PSPL) light 
curve, and in particular the short duration signatures of planets. The feasibility of such a monitoring 



campaign was demonstrated in the 1995 pilot season of PLANET (Albrow et al. 1998), during which we 
achieved ~ 2 hour sampling and few percent photometry on several concurrent bulge microlensing events. 

The MPS collaboration used observations of the high-magnification event MACHO 98-BLG-35 to rule 
out Jovian companions to the primary microlens for a large range of separations ( Rhie et al. 1999b[ ). We 
performed a similar study of OGLE-1998-BUL-14 ( Albrow et al. 2000b[ ), demonstrating that companions 
with mass > 10 Mj^p were ruled out for separations 1 — 7 AU. Our detection efficiency for this event was 
~ 60% for a companion with the mass and separation of Jupiter, thereby demonstrating that a combined 
analysis of many events of similar quality would place interesting constraints on Jovian analogs. A similar 
analysis was performed for events OGLE-1900-BUL-12 and MACHO 99-LMC-2 by the MOA collaboration 



( [Bond et al. 200l[) 



□ 



Bennett et al. (1999) claimed to detect a planet orbiting a binary microlens MACHO 97-BLG-41. As 
we discuss in §0, we exclude binaries with mass ratios q > 10^^ from our search because of the difficulty of 
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modeling binaries and therefore of making an unambiguous detection of planetary perturbations amongst 



the wealth of other perturbations that can occur in these systems. Indeed, Albrow et al. (2000a ) found that 
all available data for this event were explained by a rotating binary without a planet. 

Rhie et al. (1999b) claimed "intriguing evidence" for a planet with mass ratio 4 x 10^^ < 9 < 2 x lO""* 
in event MACHO 98-BLG-35. This perturbation had a reduced Ax^ ^ 21, far below our threshold of 60. 
As can be seen from Figure 0, our data set contains many perturbations with Ax^ < 50. As we show in 



§6.3, based on studies of constant stars, we find that systematic and statistical noise can easily give rise to 
deviations in our data with Ax^ <, 60. 



Bond et al. (2001) reanalyzed all available data for MACHO 98-BLG-35 including the then unpublished 
PLANET data that are now presented here. They found fits for 1-3 planets all with masses q < 3 x 10^^, 
with Ax^ — 60. This mass range is below our search window, primarily because our sensitivity to it is 
quite low (see §|^). In our view, planetary detections in this mass range should be held to a very rigorous 
standard, a standard not met by Ax^ = 60 which would be just at our threshold. 



Thus, none of these claimed detections ( Bennett et al. 199£ ; Rhie et al. 1999b ; Bond et al. 2001 ) would 



have survived our selection criteria even if they had been in our data. Therefore, they pose no conflict with 
the fact that we detect no planets among 43 microlensing events, and are not in conflict with the upper 
limits we place on the abundance of planets among bulge stars. 

Despite the excellent prospects for detecting planets with microlensing, and after more than five 
years of intensive monitoring of microlensing events, no unambiguous detections of Jupiter-mass lensing 
companions have been made. These null results broadly imply that such planetary companions must not 
be very common. In the remainder of this paper we quantify this conclusion by analyzing five years of 
PLANET photometry of microlensing events toward the bulge for the presence of planets orbiting the 
primary microlenses. We use strict event selection criteria to construct a well defined subsample of events. 



Employing analysis techniques presented in Gaudi fc Sackett (2000 ) and applied in Albrow et al. (2000b ) 



we search for the signals of planets in these events. We find no planetary microlensing signals. Using this 
null result, and taking into account the detection efficiencies to planetary companions for each event, we 
derive a statistical upper limit to the fraction of primary microlenses with a companion. Since most of the 
events in our sample are likely due to normal stars in the Galactic bulge, we therefore place limits on the 
fraction of stars in the bulge with planets. 

We describe our observations, data reduction and post-processing in In §|^, we describe and 
categorize our event sample. We define and apply our event selection criteria in this section also includes 
a description of how our events are fitted with a PSPL model. We summarize the characteristics of our final 
sample of events in In §|^, we describe our algorithm for searching for planetary perturbations (§^^) as 
well as various nuances in its implementation (§§ ^.2.l[ - ^.2.5| ). We describe our detections (or lack thereof) 



in §6.3 and our detection efficiencies in §3.4. Our method of correcting for finite source effects is discussed 
in §^, and we derive our upper limits in §p[ We interpret our results in §^ compare our results with other 
constraints on extrasolar planets in § p^ , and conclude in § p4] . Appendix |^ lists our excluded anomalous 
events, and Appendix ^ discusses parallax contamination. 

This paper is quite long, and some of the discussion is technical and not of interest to all readers. Those 
who want simply to understand the basic reasons why we conclude there are no planets and understand our 
resulting upper limits on companions should read §^ and §§p[]ri|. Those who want only the upper limits 
and their implications should read §§ | lO| and pi] , especially focusing on Figures ^ and |l^. A brief summary 



of this work is given in Albrow et al. (2001b) 
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2. Observations, Data Reduction, and Post-Processing 

Details of the observations, detectors, telescopes, and primary data reduction will be presented 



elsewhere (Albrow et al. 2001d). Here we will summarize the essential aspects of the observations and 



primary data reduction, and discuss only our post-processing in detail. 

The photometry of the microlensing events presented and analyzed here was taken over five bulge 
seasons starting from June of 1995 and ending in December 1999, with a few scattered baseline points taken 
in early in 2000. These data were taken with six different telescopes: the CTIO 0.9m, Yale-CTIO Im, and 
Dutch/ESO 0.91m in Chile, the SAAO Im in South Africa, the Perth 0.6m near Perth, Australia, and the 
Canopus Im in Tasmania. Measurements were taken in the broadband filters Vj and Ic using a total of 11 
different CCD detectors. 

The data are reduced as follows. Images are taken and flat-fielded in the usual way; these images are 
then photometered using the DoPHOT package ( Bchechter, Mateo, fc Saha 1993| ). A high-quality image is 



chosen for each field, which is then used to find all the objects on the frame. From this "template" image, 
geometrical transformations are found for all the other frames. Fixed-position photometry is then performed 
on all the objects in all the frames. The time-series photometry of all the objects found on the original 
template image is then archived using specialized software designed specifically for this task. This software 
enables photometry relative to an arbitrarily chosen set of reference stars. We treat each light curve for 
each site, detector and filter as independent. The number of independent light curves for each event ranges 
from one to twelve. For the majority of the events, the t^-band data are reduced using the source positions 
identified with the /-band template image, since, in general, the signal-to-noise is considerably higher in 
/-band and more objects are detected. This improves the subsequent photometry relative to what can be 
achieved using a y-band template. 

Once the photometry of all objects in the microlensing target fields are archived, we perform various 
post-reduction procedures to optimize the data quality. The light curves of the microlensing source stars are 
extracted using reference stars chosen in a uniform manner. Four to 10 reference stars are chosen that are 
close to the microlensing source star (typically within 30") and exhibit no detectable brightness variations. 
We require that the ratio of the mean DoPHOT-reported error in the measurements of each reference 
star to standard deviation all of the measurements of the star is approximately unity, with no significant 
systematic trend over the entire set of observations. Generally, the mean DoPHOT-reported error in a 
single measurement of a reference star is 0.01 mag. Reference stars are selected for each independent light 
curve, although typically the set of reference stars is similar for all observations of a particular event. Only 
those points on the microlensing event light curve with DoPHOT typesQ 11 or 13, and DoPHOT-reported 
errors < 0.4 mags are kept. Further data points are rejected based on unreliable reference star photometry 
as follows. For each reference star, the error-weighted mean is determined and the point that deviates 
most (> 3ct) from the mean is removed. The errors of the remaining points are scaled to force the 
per degree of freedom (d.o.f.) for the reference star light curve to unity. The error-weighted mean is then 
recomputed, and the entire process repeated until no > 3(t outliers remain. The outliers are reintroduced 
with error scalings determined from their parent light curves. Then, for each data point in the microlensing 
light curve, the of all the reference stars are summed. If this is larger than four times the number of 
reference stars, the data point is discarded. After this procedure, individual light curves are then examined. 



^DoPHOT types rate the quality of the photometry. DoPHOT type 11 indicates an object consistent with a point source 
star, whereas DoPHOT type 13 indicates a blend of two close stars. From our experience, all other DoPHOT types often 
provide unreliable or suspect photometry. 
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and light curves for which the microlensing target was too faint to be detected on the template image 
were eliminated. In addition, individual light curves with less than 10 points are eliminated. Since at least 
three parameters are needed to fit each light curve (see sQ), light curves with fewer than 10 points contain 
very little information. Finally, a small number (<, 10 over the entire dataset) of individual data points 
were removed by hand. These data points were clearly highly discrepant with other photometry taken 
nearly simultaneously, and were typically taken under extreme seeing and/or background conditions, or 
had obvious cosmic ray strikes near the microlensing target. Since there are only a handful of such points, 
their removal has a negligible effect on the overall sensitivity. Furthermore, these points cannot plausibly 
be produced by a real planetary signal, but would lead to spurious detections if not removed. 



3. General Considerations 

During the 1995-1999 seasons, PLANET rehed on alerts from three survey teams, EROS (1998-99), 
MACHO (1995-99), and OGLE (1995; 98-99). During these five years, several hundred events were alerted 
by the three collaborations combined. Often, there are too many to follow at one time, and PLANET must 
decide real-time which alerts to follow and which to ignore. Since the event parameters are typically poorly 
known at the time of the alert, and survey team data are sometimes unavailable, it is impossible to set forth 
a set of rigid guidelines for alert selection. The entire process is necessarily organic: decisions are made 
primarily by one (but not always the same) member of the collaboration, and secondarily by the observers 
at the telescopes, and are based on considerations such as the predicted maximum magnification and time 
scale of the event, the brightness and crowding of the source, and the number and quality of other events 
currently being followed. Our final compilation of events does not therefore represent a well defined sample. 
Some selection effects are present both in the sample of events alerted by the survey teams and the sample 
of events we choose to follow. Although these selection effects could in principle bias our conclusions, in 
practice their effects are probably quite minor, since the reasons that an event was or was not alerted 
and/or monitored (i.e. crowding conditions and/or brightness of the source, number of concurrent events, 
maximum magnification) are not related to the presence or absence of a planetary signal in the light curve. 
The one exception to this is the microlensing time scale, which as we show in is typically twice as long 
in our sample as in the parent population of microlensing events. One might imagine that, since our sample 
is biased toward longer time scale events, we are probing higher mass lenses. In fact, as we show in it 
is likely that we are primarily selecting slower, rather than more massive, lenses. Thus the bias toward 
more massive primaries is small. This is not necessarily a bias, per se, as long as we take care to specify 
the population of primary lenses around which we are searching for planets. Thus, provided that any a 
posteriori cuts we make are also not related to the presence or absence of planetary anomalies in the light 
curves, our sample should be relatively unbiased. 

We would like to define a sample of events in which we can search for and reliably identify planetary 
companions to the primary lenses. The events in this sample must have sufficient data quality and quantity 
that the nature of the underlying lensing system can be determined. Also, our method of searching for 
planetary perturbations is not easily adapted to light curves arising from non-planetary anomalies, such 
as those arising from parallax or equal mass binaries. Therefore, such events must be discarded. The 
remaining events represent the well-defined sample, which can then be search for planetary companions. 
In the next section, we describe our specific selection criteria designed to eliminate these two categories of 
events and the implementation of these criteria used to define our sample. However, for the most part, 
our events could be placed cleanly into these categories by eye, without the need of detailed modeling or 
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analysis. Examination of our full sample of light curves reveals that the events generally fall into three 
heuristic categories: 

(1) Poor-quality events. 

(2) High-quality events which are obviously deviant from the PSPL form for a large fraction of the data 

span, or are deviant from the PSPL form in a manner that is unlikely to be planetary. 

(3) High-quality events which follow the PSPL form, with no obvious departures from the PSPL form. 

(4) High-quality events which exhibit a short-duration deviation superimposed on an otherwise normal 

PSPL light curve. 



Events in the first category are the most plentiful: they consist of events with either a very small number 
of points 20), poor photometric precision, and/or incomplete light curve coverage. Events in the second 
category are those with high-quality data, in terms of photometric precision, coverage, and sampling. They 
typically consist of anomalies recognized real-time, and are comprised of both events that deviate from the 
PSPL form in a way not associated with binary lensing (i.e. finite source effects, parallax, and binary source 
events), and events arising from roughly equal- mass (mass ratio ~ 0.1) binary lenses. Events in the third 
category are high-quality, apparently normal events that follow the PSPL form without obvious deviations. 
Events in the last category are planetary candidates. 

The first two categories correspond to events that should be removed from the sample; events in 
the last two categories make up the final event sample, and should be analyzed in detail for planetary 
companions. Of course, some cases are more subtle, and the interpretation of the event is not so clear. In 
general, however, other deviations from the PSPL form are easily distinguishable from planetary deviations, 
with two caveats. First, there is no clear division between "roughly equal mass ratio" and "small mass 
ratio" binary lenses: if the mass ratio distribution of binary lenses were, e.g, uniform between q = 10~^ 
and 9=1, one would expect grossly deviant light curves, light curves with short-duration deviations, and 
everything in between. In practice, however, this does not appear to be the case, as we discuss below. 
Second, there exists a class of binary-source events that can mimic the short-duration deviations caused by 
planetary companions ( Gaudi 199§| ). Detections of short-duration anomalies must therefore be scrutinized 



for this possibility. 

All of the 126 Galactic bulgeQ microlensing events for which PLANET has acquired data during the 
1995-1999 seasons are listed in Table |l|. A cursory inspection of these events reveals that ~ 40% clearly 
belong in category (1), ~ 11% clearly belong in category (2), and ~ 25% clearly belong in category (3). 
The remaining ~ 24% are marginal events that could be placed in either category (1) or (3). However, 
no events clearly belong to the last category, i.e., there are no events that have anomalies that are clearly 
consistent with a low mass-ratio companion. Since we do not see a continuous distribution in the time 
scale of deviations with respect to the parent light curve time scale, this implies that either the mass 
ratio distribution is not uniformly distributed between equal mass and small mass ratios or our detection 



efficiency to companions drops precipitously for smaller mass ratios. In fact, as we show in \QA^ our 
efficiencies are substantial for mass ratios ~ 10^"^, implying that massive planetary companions are probably 
not typical. For the remainder of the paper, we will use strict event selection criteria and sophisticated 
methods of analysis to justify and quantify this statement. 



■^We exclude events toward the Magellanic Clouds. 
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4. Event Selection 



The goal of our selection criteria is to provide a clean sample of events for which we can reliably search 
for planetary deviations and robustly quantify the detection efficiency of companions. Such criteria are also 
necessary so that future samples of events (and possibly future detections) can be analyzed in a similar 
manner, and thus combined with the results presented here. Our selection criteria roughly correspond 
to the categorization presented in Note that any arbitrary rejection criterion is valid, as long as the 
criterion is not related the presence or absence of a planetary signal in the light curve. 

We first list our adopted rejection criteria, and then describe the criteria, our reasons for adopting 
them, and the procedure to implement them. The three rejection criteria are: 

(1) Non-planetary anomalies (including parallax, finite source, binary sources, and binaries of mass ratio 

> 0.01). 

(2) Events for which no individual light curve has 20 points or more. 

(3) Events for which the fractional uncertainty in the fitted impact impact parameter, wq, is > 50%. 



The original sample of 126 events along with an indication of which events were cut and why is tabulated 
in Table |]. The first criterion eliminates 19 events, the second 32 events, and the third 32 events, for a final 
sample of 43 events. 

As stated previously, criterion (1) is necessary because we do not have an algorithm that can 
systematically search for planetary companions in the presence of such anomalies. We are confident that the 
anomalies in the events that we have rejected by criterion (1) are, in fact, non-planetary in origin, based on 
our own analyses, analyses in the published literature, and a variety of secondary indicators. Descriptions 
of each of these events and the reasons why we believe the anomaly to be non-planetary in origin are given 
in Appendix |a| 

We fit the observed flux Fi of observatory/band / and time tk to the microlensing-event model, 

Fi{tk) = Fs,iA{tk) + Fb,i + m[0{tk) - ^0.;] (6) 

where A(tk) is the magnification at time tk] ^s,i and i<B,; are the source and blend fluxes for light curve 
I. The last term is introduced to account for the correlation of the flux with seeing that we observe in 



almost all of our photometry (see Albrow et al. 2000b). Here rji is the slope of the seeing correlation, 9{tk) 



is the full width at half maximum (FWHM) of the point spread function (PSF) at time tk, and Oqj is the 
error- weighted mean FWHM of all observations in light curve /. For a single lens, the magnification is given 



by ( [Refsdal 1964| ; [Paczynski 19861 ) 



u{t)^^ 

where t is the "normalized time," 



Ao[u{t)]= : u'{t)=r^+ul (7) 



t-tn 



(8) 



Here to is the time of maximum magnification, is the characteristic time scale of the event, and uq is 
the minimum angular separation (impact parameter) between the lens and source in units of 9e- A single 
lens fit to a multi-site, multi-band light curve is thus a function of 3 -I- 3A^; parameters: ^e, wo, ^o, and one 
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source flux fs.ij blend flux -Fb,;, and seeing correlation slope 77; for each of Ni independent light curves. For 
a binary lens, three additional parameters are required: the mass ratio of the two components, q, the binary 
separation d in units of 9^, and the angle of the source trajectory with respect to the binary axis, a. Thus 
for an event to contain more information than the number of free parameters, at least one observatory must 
have at least 9+1=10 data points. In order for the fit to be well-constrained, considerably more data points 
than fit parameters are needed. We therefore impose criterion (2): if no independent light curve has at least 
20 data points, the event is rejected. The number 20 is somewhat arbitrary, however the exact choice has 
little effect on our conclusions: a natural break exists such that the majority of events are well above this 
criterion, and those few events that are near the cut have little sensitivity to planetary perturbations. 

All events that pass criterion (2) are fit to a PSPL model [eqs. and (Q)]. At this stage, we also 
incorporate MACHO and/or OGLE data into the fit, when availablen. To fit the PSPL model, we combine 



the downhill-simplex minimization routine AMOEBA ( Press et al. 199^ ) with linear least-squares fitting. 
Each trial combination of the parameters (iEi^Ei^io) immediately yields a prediction for AQ{t) [eqs. (|^) 
and (H)]. The flux is then just a linear combination of Fg^i, Fb^i and rji [eq. (^)]. The best fit parameters 
a-i = (Fs,i,^b,i,»71)-Ps,2,Fb,2,??2, •■•) can then be found by forming, 

da. da, ' ' da. ' 

where the index k refers to a single observation, the sum is over all observations, and ak is the photometric 
error in the observed flux F^. The parameter combination that minimizes is then, 

j 

Occasionally, the values of i^B,; obtained from this procedure are negative. If i^e,; is negative by more than 
its uncertainty, we apply a constraint to to force Fbj = 0. We then use AMOEBA to find the values of 
{tE,tE,uo) that minimize x^- Note that since neither MACHO nor OGLE report seeing values, we do not 
correct their data for seeing correlations. 



We know from experience ( [Albrow et al. 1998| , pOOOq ) that DoPHOT-reported photometric errors are 



typically underestimated by a factor of ^ 1.5. Naively adopting the DoPHOT-reported errors would thus 
lead one to underestimate the uncertainty on fitted parameters, and overestimate the significance of any 
detection. However, simply scaling all errors by a factor to force x^/d.o.f. to unity is also not appropriate, 
as we find that our photometry usually contains significantly more large (> 3cr) outliers than would be 



expected from a Gaussian distribution (Albrow et al. 20001, 2001a). Furthermore, independent light curves 



from different sites, detectors, and filters typically have different error scalings. Therefore we adopt the 
following iterative procedure, similar to that used by Albrow et al. (2000b| ). We first fit the entire dataset 



for a given event to a PSPL model in the manner explained above. We find the largest > Str outlier, and 
reject it. We then renormalize the errors on each individual light curve to force x^/d.o.f. to be equal 
to unity for that light curve. Next, we refit the PSPL model, find the largest > 3cr outlier, etc. This 
process is repeated until no > Su outliers are found. The outliers are then reintroduced, with error scalings 
appropriate to their parent light curve. We typically find 3 to 6 outliers > 3ct in the PLANET data and 
OGLE data, and a larger number for MACHO data (which contain significantly more data points). The 



^MACHO data are available for those events alerted by MACHO in 1999, along with a few events that were originally alerted 
by OGLE in 1999. OGLE data is available for events alerted by OGLE in 1998-99, along with a few events that were originally 
alerted by MACHO during these years. 
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median error scaling for PLANET data is 1.4, with 90% of our data having scahngs between 0.8 and 2.8. 
The errors reported by OGLE are typically quite close to correct (scalings of ^ 1.1), while MACHO errors 
are typically overestimated (scalings of ~ 0.8). 

Once the best-fit PSPL model is found, we determine the uncertainties on the model 
parameters by forming Cy as in equation (^), except that now the parameters are 
(to, ^E, uq, -Fs,i5 -Fb,i, ^71, ^s,2, ^B,2, 7?2, •••): i-^-j we have included t-E, to, and uq- The uncertainty in 
parameter is then simply Sai — {cuY^'^. Note that we include the outliers to determine the uncertainties. 



As discussed by Griest fc Safizadeh (1998 ) the sensitivity of a light curve to planetary companions is 
strongly dependent on the path of the source trajectory in the Einstein ring, such that trajectories that 
pass closest to the primary lens, i.e. events with small uq, will have larger sensitivity than events with 
larger uq. Thus, in order to accurately determine the detection efficiency to a given binary lens, the source 
path in the Einstein ring, u(t), must be well-constrained; poor knowledge of u{t) translates directly into 



poor knowledge of the sensitivity of the event to planets ( Gaudi fc Sackett 200C ). The values of u{t) for 
a given dataset are determined from the mapping between flux and magnification, which depends on the 
source and blend fluxes, and the mapping between the magnification and time, which depends on uq, t^, 
and to. In blended PSPL fits, all these parameters are highly correlated. Thus, a large uncertainty in uq 
implies a large uncertainty in other parameters. Thus the uncertainty in uq in a PSPL fit can be used as an 
indication of the uncertainty in u(t), and thus the uncertainty in the detection efhciency. Furthermore, for a 
planetary perturbation, the projected separation c? is a function of the observables (^o.p — to) /t-E, where to.p 



is the time of the planetary perturbation, while the mass ratio is q ^ tp/^E ( Gould fc Loeb 1992 ; Gaudi fc 
Gould 1997 ), where tp is the duration of the perturbation. Therefore the detection of a planet in an event 



with poorly constrained t^iuo) would be highly ambiguous, as the neither the projected separation d nor 
the mass ratio q would be well-constrained. We therefore impose a cut based on the fractional uncertainty 
in the fitted value of uq. 

Figure |^ shows the fractional uncertainty Suq / uq in the impact parameter versus uq for all events that 
passed selection criteria (1) and (2). Examination of the distribution of fractional uncertainty in uq for 
these events reveals a large clump of events with small fractional uncertainty; many scattered, smoothly 
distributed events with larger uncertainties, and a natural break in the distribution at Suq/uq ~ 50%. We 
therefore adopt Suq/uq = 50% for our final event cut. The exact choice for the cut on Suq/uq has little 



effect on our conclusions; as we discuss in § |6.4 events with Suq/uq ^ 30% typically have low detection 
efhciencies. Four classes of events have poorly-constrained uq. These are events: for which the data cover 
only one (usually the falling) side of the event; for which no baseline information is available; that are 
highly blended; with an intrinsically low maximum magnification. Thus by imposing a cut on Suq/uq, we 
eliminate all low magnification events; the event with largest impact parameter in our final sample has 
Mo = 0.61. Note that the majority of events that fail the cut on Suq/uq fall into the first two classes, which 
emphasizes the need for coverage of the peak and baseline information. In particular, without MACHO and 
OGLE data, many more events would not have passed this last cut, and our final sample would have been 
considerably smaller. 

After imposing cuts 1 (non-planetary anomalies), 2 (data quantity), and 3 (uncertainty in the impact 
parameter), we are left with a sample of 43 events. The light curves for these events are shown in Figure ^. 
In order to display all independent light curves (which in general have different Fg, Fq, and 77), we plot 
the magnification, which is obtained by solving equation (^) for Ao{t). Rather than show the magnification 
as a function of true time, we show the magnification as a function of normalized time r [eq. (j^)]. When 
plotted this way, perturbations arising from a given q would have the same duration on all plots [eq. (j^]. 
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Thus the sensitivity of different Hght curves to companions can be compared directly. In the next section, 
we describe the properties of these events, paying particular attention to those properties relevant to the 
detection of planetary anomalies. 



5. Event Characteristics 

The parameters to, ^e, and uq and their respective la uncertainties for the final event sample are 
tabulated in Table 0, along with the percent uncertainty in uq. The sensitivity of an event to planetary 



companions depends strongly on uq (Gould fc Loeb 1992; Griest fc Safizadeh 1998; Gaudi fc Sackett 2000) 



and thus the exact distribution of uq influences the overall sensitivity of any set of light curves. The time 
scale is important in that the population of lenses we are probing is determined from the distribution of t^- 
In addition, we use in §0 to estimate the effect of finite sources on planetary detection efficiencies and 
therefore the effect on our final conclusions. For the current analysis, the parameter to is of no interest. 

In Figure ^, we plot uq against t-E for our event sample, revealing no obvious correlation between the 
two. This lack of correlation between and uq implies that the lenses that give rise to the events with the 
most sensitivity to planets (i.e., those with small uq) comprise a sample that is unbiased with respect to the 
entire sample of lenses. Given this, we can then inspect the distributions of uq and ^e independently. 

Both the differential and cumulative distributions of ^e are shown in Figure ^. The median time scale 
of our events is ~ 40 days, about a factor of two higher than the median time scale for events found by 
the MACHO and OGLE teams toward the Galactic bulge ( [Alcock et al. 1997a| ; [Udalski et al. 2000| ). This 



is almost certainly a selection effect caused by the fact that longer time scale events are more likely to be 
alerted before peak magnification, and thus are more likely to be chosen by us as targets for follow-up 
photometry. This is compounded by the fact that, for short time scale events, we are less likely to get good 
coverage of the peak, even if they are alerted pre-peak. Events with poor or no peak coverage will often fail 
our selection criterion of < 50% fractional uncertainty in uq. In principle, this deficiency could be partially 
alleviated by including MACHO and/or OGLE data. However, in practice, we often stop observing the 
event altogether if we do not get good peak coverage. As we discuss in the primary effect of this selection 
is a bias toward slower lenses. 

We also show in Figure ||the differential and cumulative distributions of uq. The median uq is ~ 0.2, 
and the fraction of high-magnification (^max > 10) events is ~ 30%. As it is a purely random quantity, 
the intrinsic distribution of uq should be uniform. The observed distribution of uq, however, is clearly not 
uniform. This is due to a combination of various selection effects. First, faint events are more likely to be 



detected (and hence alerted) by the survey teams if they have a larger maximum magnification ( Alcock 



et al. 1997a| ; [Udalski et al. 20001 ). Since there are more faint stars than bright stars, this results in a bias 



toward smaller impact parameters with respect to a uniform distribution. Second, since events with smaller 
impact parameters are also more sensitive to planets, we preferentially monitor high-magnification events. 
This bias does not affect our conclusions, since the value uq is unrelated to the presence or absence of a 



planetary companion. However, as emphasized by Gaudi fc Sackett (200C ), it does imply that in order to 



determine accurately the overall sensitivity of an ensemble of light curves to planetary companions, the 
actual distribution of observed uq must be used. 

Since one of the primary goals of PLANET is to obtain very dense sampling of microlensing events, 
it is interesting to examine how well this goal has been achieved. In Figure ^ we show the distribution 
of sampling intervals, that is, the time between successive exposures of a given event. Three peaks are 
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evident. The first at ~ 6 minutes is our typical J-band exposure time of 5 minutes plus 1 minute of 
overhead time; this peak is dominated by events that are followed continuously and also pairs of I-V data 
points. The second peak at ~ 1.5 hours represents our fiducial sampling interval. The third peak at 1 day 
arises primarily from sampling of the wings and baselines of light curves. The median sampling interval is 
~ 1.5 hours, with 90% of all data taken between 5 minutes and 1 day of one another for a given event. What 
is of particular relevance to the detection of planets is the sampling interval in units of t-^, which is shown in 
the lower panel of Figure ^. Assuming that at least 10 data points are needed on a planetary perturbation 
for detection, the sampling interval needed to detect a companion of mass ratio q is approximately, 



Ai = 3x lO-^tE^^. (11) 

Using this formula and comparing to Figure g we find that (80%, 65%, 45%, 25%) of our data have sufficient 
sampling to detect companions of mass ratio (lO^'^, 10^'^, 10^**, 10~^). Thus we expect the majority of our 
data to have sufficient sampling to detect companions with mass ratios > lO^"'. This is not an accident, 
since PLANET observations are planned to have sensitivity to Jovian mass planets orbiting main sequence 
stars ( lAlbrow et al. 199S| ). 



The sensitivity of a given light curve to planetary companions is primarily determined by three factors: 
photometric errors, temporal sampling, and impact parameter. In Figure |^, we plot the median photometric 
error, CTmcd, versus the median sampling interval, A^med for all events; high-magnification (Amax > 10) events 
are indicated. These are also tabulated in Table ||. High magnification events that occupy the lower left 
quadrant of Figure ^ will have the highest sensitivity to planetary companions. Of the 13 high-magnification 
events, all have sufficiently small median sampling intervals to detect q — 10^^ companions; we therefore 
expect our sensitivity to such to companions to be quite high. Two high-magnification events have sufficient 
sampling rates to detect companions with q = 10"^; however, for companions as small as this, excellent 
photometry (< 2%) along with excellent sampling is required to obtain significant efficiency for detection 
( Bennett fc Rhie 1996 ). No events satisfy both of these requirements ((Tmed < 2% and Atmcd/^E < lO^'^ '^). 



We therefore restrict our attention to g > 10 . 

Considering the large number of high-magnification events, and the dense sampling and precise 
photometry, our sample should be quite sensitive to planetary companions, especially those with q > 10~^. 
This fact, combined with the fact that no planetary- like perturbations are clearly evident in the light curves, 
is an indication that such planetary companions are probably not common. In the following sections, we 
strengthen and quantify this statement. 



6. Search for Detections and Calculation of Detection EfRciencies 

Although a cursory inspection of Figure |^ reveals no obvious candidate planetary perturbations, 
such perturbations could be quite subtle, and thus missed by eye. Furthermore, the significance of the 
lack of planetary perturbations must be quantified. Specifically, the frequency with which planetary 
companions of given d, q could be detected in individual light curves, the detection efficiency, must be 
determined. We simultaneously search for planetary signatures in and determine the detection efficiency of 
individual events using the method suggested by Gaudi fc Sackett (200(]| ) and applied to microlensing event 



OGLE-1998-BUL-14 by [Albrow et al. (2000bD . We briefly review the algorithm here, but point the reader 



to these two papers for a more thorough discussion of the method and its application. 



- 14 - 



6.1. Algorithm 

Of the 6 + 3iV; parameters (see in a point source binary microlensing fit, 3 + 3-/V; have analogs in 
the PSPL fit: tg, uq, <o a-nd one Fg F^.i, and ?// for each of Ni independent Hght curves. The parameters 
tE, uo and to have different meanings in the binary- lens model than in the PSPL model, and depend on the 
choice of the origin of the binary-lens and the reference mass. For small mass-ratio binaries, however, if 
one chooses the origin to be the location of the primary lens, and normalizes to the mass of the primary, 
then the values of these parameters will be quite similar in a binary-lens and single-lens fit to a light 
curve. Three parameters are not included in the PSPL fit: the mass ratio q, the projected separation d, 
and the angle a of the source relative to the binary- lens axis. While q and d are related to the physical 
nature of the planet-star system, the angle a is a nuisance parameter which is of no physical interest. It 
is a random geometric parameter and therefore uniformly distributed. However, the value of a does have 
a significant effect on the amplitude and duration of the planetary perturbation. Thus, some values of a 
lead to detectable perturbations to the PSPL model, while others do not. Marginalization over a for a 
given binary lens specified by {q,d) therefore determines the geometric detection efficiency ei(d,q) for event 
i and such a binary system. Repeating this process for all (d, q) pairs of interest yields the efficiency for 
all systems. This is the basis of the method of determining the detection efficiency for individual events 
suggested by Gaudi fc Sackett (2000|) . 



Operationally, the procedure to search systematically for planetary signatures and determine for 
each event is as follows: 

(1) Fit event i to the PSPL model, obtaining x|spl (§@)- 

(2) Holding d and q fixed, find the binary-lens model that best fits light curve i for source trajectory a, 

leaving the 3 -I- 37V; parameters (ts, uo,to, [Fg , i^B , '7] TV/ ) as free parameters. This yields x^id, q, «)• 

(3) Repeat step (2) for all source trajectories < a < 27r. 

(4) Evaluate the difference in between the binary and PSPL fits: Ax^(d, g, a) = x^(d, g, a) — XpgpL. 

Compare this to some threshold value Axthrcsh- 

(a) If Ax^{d,q,a) < — Axthj-cgh' ^1^™ tentatively conclude we have a detected a planet with 
parameters d, q, and a. 

(b) If Ax^(d, <z, a) > Axthresh then the geometry {d,q,a) is excluded. 

(5) The detection efficiency ei{d, q) of event i for the assumed separation and mass ratio is then 



<d, q)^^ j^J da e[Ax2(d, (?, a) - AxLesh], (12) 

where Q[x\ is a step function. 

(6) Repeat steps (2)-(5) for a grid of {d,q) values. This gives the detection efficiency ei{d,q) for event i as 
a function of d and q, and also yields all binary-lens parameters (c?, q, a) that give rise to significantly 
better fits to the event than the PSPL model. 



(7) Repeat steps (l)-(6) for all events in the sample. 
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In step (2), we find the parameters (tE, wq, io? [^Sj ^b, '7]-^/ ) that minimize in the same way as the 
PSPL fit: we choose trial values of (tE, wo, ^o) which (along with the values of d, q, a) immediately yield the 
binary- lens magnification^ as a function of time, A^(£). This is used to find the least-squares solution for 
the other parameters, and the resultant x^- A downhill-simplex routine is then used find the combination 
of parameters (tE,wo,to) that minimize (see §||). The procedure is slightly more complicated for those 
events for which MACHO and/or OGLE data was used for the PSPL fit, as we discuss in § |6.2.2 . 

Due to the perturbative nature of the planetary companion, for the appropriate choice of the origin 
of the binary and the total mass of the system, the majority of structure of the hypersurface with 
respect to the parameters (iE, uq, io) will be very similar in the PSPL and the binary lens cases. The two 
hypersurfaces will only deviate significantly in some localized region of the (iE,MO:^o) parameter space 
where the planetary perturbation from the PSPL form is large. Consider a set of parameters (fi, (7, a) 
for which the characteristic size of such a region in (tE,wo,to) space is much smaller than the intrinsic 
uncertainty of these parameters. Since we find the binary- lens fit that minimizes x^j rather than integrating 
over the whole x^ surface, our algorithm will find best-fit parameters (tE,uo,to) for the binary-lens model 
that avoids this region without significantly increasing the x^ with respect to the single lens. Thus we 
will always underestimate the detection efficiency. The amount the detection efficiency is underestimated 
depends on how well iE,uo, and are constrained. For events with poorly-constrained parameters, the 
efficiency can be underestimated by a significant amount ( Gaudi fc Sackett 2000 ). This is illustrated in 
Figure ^ using event OGLE-1998-BUL-13 as an example. The fractional uncertainty in uq for this event 
is ~ 7%. We show the vector positions in the source plane of the data points for this event for the best-fit 
as determined from the PSPL fit, along with the ±4(T bounds on uo0 The data are more "compressed" 
in the Einstein ring for values of uq smaller than the best-fit value because ^e is anti-correlated with uo, 
and thus smaller wq implies larger ^e- For reference, we also show contours of constant fractional deviation 
from a single lens for a binary with q = 0.001 and 6=1.11. It is clear that the difference in between the 
binary-lens and single-lens fits will differ substantially between these three fits. Our algorithm will always 
choose the one that minimizes x^i thus will underestimate the efficiency. This could in principle be 
avoided by integrating over uq, to, and ^e, rather than evaluating x^ at the best-fit parameters. However, for 
the large number of binary-lens geometries we test (see § |6.2.3 ), this is not computationally feasible. These 
underestimated detection efficiencies could be a serious problem if planetary deviations were detected, as 



they would lead to an overestimate of the true number of planets. However, as we show in §6.3, we do not 



detect any planetary deviations. Thus, the underestimated efficiencies represent conservative upper limits. 



6.2. Implementation of the Algorithm 



Although the algorithm described in §6.1 is conceptually simple and appears straightforward, there are 
some subtle details that must be addressed before implementation. Specifically, in the following subsections 
we discuss photometric errors, the inclusion of MACHO /OGLE photometry, the grid size and spacing for 
the binary-lens parameters d, q, and a, the method by which the binary-lens magnification is evaluated, 
and the choice of the detection threshold Axthrcsh- 



For an explanation of how to calculate the binary-lens magnification, see Witt (199C ). 

^Note that the bounds on ug were calculated by projecting the Ax^ surface on uq, rather than by the linearized covariance 
matrix, as in Table ^. In general, the former method gives asymmetric bounds on uq due to the Fg > constraint, whereas 
the latter gives symmetric bounds by definition. 
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6.2.1. Photometric Errors 

As we discussed in the errors reported by DoPHOT are typically underestimated by a factor of 
~ 1.5; adopting such errors would both overestimate the significance of any planetary detections, and 
overestimate the detection efficiency. Furthermore, since events can have error scaling factors that differ by 
a factor of three, even the relative significances for different events would not be secure. Ideally, one would 
like to determine the magnitude of the photometric errors without reference to any model. Unfortunately, 
this is not possible in general, primarily because the error depends strongly on the local crowding conditions 
of the microlensing source object in a manner that is impossible to access a priori. Therefore, in order to 
put all events on the same footing and to arrive at the best possible estimate of the significance of planetary 
detections and detection efficiencies, we adopt the error scaling factors as determined in the PSPL fit (see 
§H). We typically find that, after scaling in this way, the error distributions are nearly Gaussian, with the 



exception of a small handful of large > 3ct outliers ([Albrow et al. 2000q ) 



If the PSPL model is truly the "correct" model, this procedure is valid, and does not bias the results 
in any way. However, if the light curve actually deviates from the PSPL model, this procedure will 
overestimate the error scaling factors, and thus underestimate the significance of the anomaly. Assuming 
that binary-lens model is correct, it is straightforward to show that the true difference in x^, which we will 
label Axg, is related to the Ax^ evaluated assuming the PSPL fit is correct by, 

where d.o.f. is the number of degrees-of-freedom of the event. Thus for an event with 300 data points and 
Ax^ — 60, using the errors determined from the PSPL fit would lead us to underestimate the "true" Axo 
by 20%. For events with a small number of d.o.f., this underestimate can formally be as large as 100%. This 
would seem to argue that the values of x^ computed in all fits (PSPL and binary) should be renormalized 
by the best- fit model (PSPL or binary). However, there are several reasons we feel this is not appropriate. 
First, for any fit, x^ is not dominated by the number of d.o.f.: instead, typically only a handful of large 
outliers contribute a significant fraction of the evaluated x^- Thus, in reality d.o.f. should be replaced by 
^binary cquatiou ([T^) , which is typically larger by 100, thus reducing the underestimate considerably. 
Furthermore, renormalizing x^ in this way would give extra weight to binary-lens models that "succeed" 
by fitting isolated large-a outliers, particularly for events with a small number of data points, where x^ is 
dominated by such outliers. The smaller the number of data points, the more difficult it is to objectively 
judge the reality of such fits. Although some of these biases could in principle be calibrated by Monte Carlo 
techniques, i.e. by inserting many artificial planetary signals into the light curves, and then repeating the 
algorithm on all of these artificial datasets, in practice the large number of fits required (see § |6.2.3 ) makes 



this computationally prohibitive. Furthermore, it is difficult to address the effects of large-tr outliers in 
this way. We will therefore adopt the conservative and simpler choice of using the errors determined with 
reference to the PSPL model in order to avoid the danger of detecting spurious planets in data with isolated 
outliers in sparse datasets. 



6.2.2. Including MACHO/OGLE Data 

As discussed in we include MACHO and/or OGLE data for some events in order to better constrain 
uq. This is necessary in order to robustly determine ti for events for which our data are poorly sampled near 
the peak or do not have baseline information. However, as we do not have access to these raw data, nor do 
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we know the details of the data reduction procedures, we have no way of independently judging the quality 
of the MACHO or OGLE photometry. Furthermore, we do not have access to the seeing values for these 
data, and hence cannot correct for the seeing correlations that can often mimic low-amplitude planetary 
deviations. Thus any planetary "signal" discovered using this photometry would be difficult to interpret, 
and the reality of the signal impossible to determine. Therefore, while we use these data to constrain the 
global parameters ^e, io and mq, we do not use these data in either the search for planetary signatures or 
the calculation of the planet detection efficiency. We accomplish these goals in the following manner. 

All information on the parameters t-^, to, uq and their covariances with other parameters is contained 
within the covariance matrix Cij and the vector di as determined from the PSPL fit with all parameters [see 
^ and eq. (^]. Therefore, we simply need to extract the information provided by the MACHO/OGLE data 
and apply it to the binary-lens fit with only PLANET data. First we calculate the covariance matrix Cij 
of the best-fit parameters — (io: ^Ei i^Oj ^S,ii ^B.ii ^ii ^S,2i ^b,2: ••■) as determined by the PSPL fit to 
all (MACHO+OGLE+PLANET) data. Note that this is identical to the procedure used in §| to calculate 
the uncertainties of ai. We then restrict Cij and Ui to the parameters Fg^, Fbj,'!]! for PLANET data. We 
call these restricted quantities c^°^ and af^^^ . We calculate the covariance matrix of the best-fit 
parameters of determined from the PSPL fit to only PLANET data, again restricting these quantities to 
the parameters -Fs,;j Fb,^ Vi- Next, we form the matrix and vector, 

^MOP ^ (c^OP) -1 ^MOP = J2 6^°Paf (14) 

j 

and similarly for bfj and df . Finally, we calculate, 

6^o = 6^oP-6P dr=dfOP-df. (15) 

The resultant matrix and vector d^*-* contain only the information on to,tE,uo and the parameters 
Fs,i, FB,i,r]i for PLANET data provided by the MACHO/OGLE data. We then use these two quantities to 
constrain the binary-lens fits using PLANET data only in the following manner. For each trial ttj,t-E,UQ, 
we compute bij and di for the quantities Fs.i, FBj,rii using only PLANET data. We add to these the 
constraints from MACHO/OGLE by forming 

b'^°^' = b,j+b^° d™"'' = d, (16) 

which are then used to find the best-fit parameters = (-Fs.ij ^b,i, ?yi, -Fs,2, -Fb.2, ^72, ••■) via equation ([l0|). 
The of the resultant fit is then evaluated. We add to this a contribution, 

where = (6^0) -i ^nd 

5a. =a.- af O, af « = ^ c^^^df (18) 

j 

The contribution Xmo ^ ^ penalty for violating the constraints from MACHO/OGLE data. The 
remainder of the fitting procedure is as before: this x^ is then used by the downhill-simplex routine 
AMOEBA ( Press et al. 1992j ) to find the parameters and that minimize x^ for the particular d, q, a 



binary-lens geometry. 
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6.2.3. Grid o J Binary- Lens Parameters 



Several factors dictate our choice of grid size and spacing in d, q, a parameter space. First, the grid 
spacing must be dense enough to avoid missing possible planetary signals and prevent sampling errors from 
dominating the uncertainty in e^. Second, the grid must cover the full range of parameter space for which 
we have significant sensitivity. Finally, the computation must be performed in a reasonable amount of time. 

We restrict our attention to lO"** < q < 10^^. The upper end of this range is dictated by the fact that 
we are primarily interested in planetary companions, and also because our procedure for finding binary-lens 
fits fails for events that are grossly deviant from the PSPL form. In fact, finding all satisfactory fits to such 



binary- lens light curves is quite difficult (see Mao fc Di Stcfano 1995; Pi Stefano fc Perna 1997; Albrow 



et al. 1999b| ). We do detect binaries well fit by g > 0.01. Incorporating such binaries into the analysis 
would entail finding all possible fits to these observed binaries and calculating the efficiency of all other 
events. Although such a study is interesting in its own right, it would be quite an undertaking, well beyond 
the scope of this paper. The lower end of the range of mass ratios we test is dictated by the fact that we 
are unlikely to have significant sensitivity below q = 10~* (§||). We sample q at equally spaced logarithmic 
intervals of 0.25. 



Numerous studies (Gould fc Loeb 1992; Di Stefano fc Mao 1996; Bennett fc Rhie 1996; Griest & 



Safizadeh 1998; Rhie et al. 1999b; Albrow et al. 2000b) have shown that the planetary detection probability 
is largest in the "lensing zone," 0.6 < d < 1.6, and is negligible for d <^ 0.1 and d ^ 10. Furthermore 
planetary perturbations exhibit a c? — > symmetry (Gaudi fc Gould 1997; Griest fc Safizadeh 199g; 



Dominik 1999b| ). Therefore, we sample d at 0.1, 0.2, 0.9, 1.0, and also the inverse of these values, for a 



fuU range of 0.1 < d < 10. 

In order to avoid missing any possible planetary signals, we choose a variable step size for a that 
depends on q. The size of the region of significant perturbation is ~ 9^^^, and thus a perturbation at the 
Einstein ring radius would cover an opening angle with respect to the center of the primary lens of ~ g^/^. 
Therefore in order to sample the perturbed region at least twice, we choose a step size of 



Aa 



2 ■ 



(19) 



For every d,q pair, we thus find the best-fit binary-lens model for a total of iirq ~ 400(^/10 ■^) 
choices of a. 



6.2.4. Magnification Maps 



With the grid size and spacing described in § |6.2.3 , we perform a total of 8.8 x 10"* binary-lens fits to 



each event, for a grand total of 3.8 x 10^ fits for all 43 events . Each fit requires at least 50 evaluations of 
the binary-lens magnification light curve to converge, for a total of more than 10® binary-lens light curve 
evaluations. Given this large number of evaluations, re-evaluating the magnification for each data point 
of each event is both prohibitive and inefficient. We therefore first create magnification maps for each of 
the d, q grid points, and interpolate between these maps to evaluate the binary-lens magnification. Maps 
are generated for source positions — 2 < a; < 2 and — 2 < y < 2 (in units of 9e). For source positions 
outside this range, we use the PSPL magnification. For a binary with q 1 and d 1, there are two 
sets of caustics. The "central caustic" is always located at the position of the primary, i.e. x — 0, y = 0. 
The "planetary caustic(s)" are separated from the primary by an amount \d^^ — d\. Therefore by only 
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evaluating the binary-lcns magnification for source positions in the ranges above, we are imphcitly assuming 
that we are not sensitive to the planetary caustics of companions with separations d <^ 0.4 and d ^ 2.4, 
although we are still sensitive to such planets via the central caustic. This assumption is essentially correct 
since the vast majority (~ 95%) of the data was taken within < of the peak. To generate the maps, 
the source position is sampled at intervals of 2 x 10~'^6'e, the typical sampling interval of our events 
We have performed numerous tests comparing fits using these maps and fits using the exact binary-lens 
magnification, and find that using the maps introduces an error of Ax^ <, 2, which is far below any of our 
thresholds Axjj^j.^j,jj. Typically, efficiencies determined using these maps are in error by <^ 1%. We have also 
inserted planetary deviations into selected light curves, and confirm that these "detections" are recovered 
when the maps are used to evaluate the magnification. 



6.2.5. Choice of Detection Threshold 

Ideally, one would like to choose the detection threshold Ax^Yircsii ^ priori, without reference to the 
results of the binary-lens fits. Specifically, one would like to be able to determine the probability P(> A^^) 
of obtaining a given A^^ or larger by chance, and then choose a probability threshold for detection, say 
P ~ 0.01. Naively, one might expect that the probability of getting a certain value of Ax^ or larger by 
chance is given by, 

/>oo 

P(> Ax') = (2^)"'/' / dx xi/2e--/2, (20) 

for the three extra binary parameters (d, q, a), assuming they are independent and have Gaussian distributed 
uncertainties. However, this formula fails for several reasons. First, most events contain large outliers that 
are not described by Gaussian statistics. Second, and more importantly, such a naive calculation fails to 
take into account the fact that many independent trial binary-lens fits to the datasets are being performed, 
thereby effectively increasing the difference in the number of degrees-of-freedom between the binary and 
single lens models. In other words, while the success of a single binary lens model is given by equation ( p0| ) 
in the limit of Gaussian errors, the success of any binary- lens model is not. Unfortunately, the effect of this 
increase in the effective number of degrees-of-freedom on the probability cannot be assessed analytically, 
and must be determined via a Monte Carlo simulation. This would entail generating many different 
realizations of synthetic events with sampling and errors drawn from the sampling and error distributions 



of each of the 43 events in our sample. The algorithm in §6J would then need to be performed on each 
of these synthetic events, in order to determine the mapping P{> Ax') for each event. Given that each 



event requires ~ 10^ binary-lens fits, this is clearly impossible. Furthermore, as we demonstrate §6.3, it 
is likely that unrecognized systematics exist in the data which give rise to temporal correlations in the 
fluxes of observed light curves. These systematics will result in false detections. The rate of such false 
detections cannot be recovered with Monte Carlo simulations of synthetic light curves unless the actual 
temporal correlations (which arc not understood) are introduced in these light curves. We therefore use the 

2 

thresh ' 



distribution of Ax' from the actual events to choose Ax^. ,^ described in the next section 



6.3. Detection Threshold and Candidate Detections 

We have applied the algorithm presented in |6.l| for all 43 events in our final sample. For each event, 
we find the absolute minimum Ay^- from this procedure. The distribution of these Ay^,- is shown in 
Figure H. If all the events harbored planets, we would expect a continuous distribution in Ax' extending 
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to very large negative values. If some fraction of events harbored planets, then we would expect a large 
"clump" of small A^^^jj^ obtained from single events through statistical fluctuations, and then a few 
scattered instances of large A^j^^^jj^ from those events with companions. In fact most of the events have 
^Xmin ^ -60, with only two events, MACHO 99-BLG-18 and OGLE-1999-BUL-36, having Axf^i^ < -60. 
We therefore interpret the binary- lens fits with Ax^j^ > — 60 to be arising from statistical fiuctuations or 
unrecognized low-level systematics, and choose Ax^jjj.^gj^ = 60 as a reasonable threshold for detections. 

To establish the plausibility of our choice of Ax^j^j,^gj^, we perform a simplistic Monte Carlo simulation. 
For one observatory and filter, we extract 1000 light curves of stars in the field of a typical microlensing 
event. These stars span a large range of brightness and local crowding conditions. The overwhelming 
majority of these stars have constant brightness, although a handful are almost certainly variables. We 
reduce and post-process these light curves in the same manner as the microlensing events (§H), using a 
constant flux model with seeing correlation correction to rescale the errors. Outliers (> 3cr) are included, 
but not used to determine the error scaling. We then fit each of these light curves to the model designed to 
mimic the deviation induced by a planetary companion: 

F{tk)^ Fs[l + 5oexpi~Ti)]+T^[eitk)-9o], Tk = (tk - to)/tp. (21) 

This model has a deviation from constant fiux with a maximum amplitude of So at a time tg, and a 
characteristic duration tp. We vary Sq in 80 steps 6o = —20% to 20%, to in 30 steps between the minimum 
and maximum time of observations, and tp in 30 logarithmic steps between 10^^ and 10^** of the total 
duration of the observations, for a total of 7.2 x 10^ trial combinations. This is similar to the number of 
binary-lens fits performed for each microlensing event. For each 6o,to, and tp, we find the best-fit values of 
Fs and 77, and calculate x^- This is repeated for all sampled values of {So,to,tp) and the minimum A^^j^^ 
between the best fit to the model in equation ( pT| ) and the constant flux model determined for each of the 
1000 fight curves. In Figure ^ we show the resulting distribution of Ax^j^^, normalized to 43 events. The 
similarity to the distribution of Ax^in of the microlensing events is remarkable. We conclude that it is 
quite likely that the binary- lens flts with Ax^in > — 60 arise from statistical fluctuations or unrecognized 
low-level systematics, and that our choice of Axt^rcsh reasonable. 

Based on this choice of Axthrcsh — ^0, we tentatively conclude that we have detected anomalies 
consistent with planetary deviations in events MACHO 99-BLG-18 and OGLE-1999-BUL-36. We have 
examined both events individually, and find other, more likely, explanations for their anomalous behavior 
which wc now describe in some detail. 

The light curve for OGLE-1999-BUL-36 shows an overall asymmetry will respect to the time of 
maximum magnification. This asymmetry is well fit by the distortion to the overall light curve created 
by a planetary companion to the primary lens with q = 0.003. However, such a distortion requires a 
special geometry, specifically a ~ or 180°, i.e. a source trajectory nearly parallel to the planet-star 
axis. All other values of a produce either no asymmetry or a planetary "bump." Asymmetries like that 



of OGLE-1999-BUL-36 are a generic feature of low-amplitude parallax effects (Gould, Miralda-Escude 



Bahcall 1994); indeed the event is equally well-fit by a parallax model. Typically, parallax effects are only 
significant in long time scale events (iE ~ 100 days), and thus it would seem unlikely that, for typical lens 
masses and distances, such effects should be detectable in the light curve of OGLE-1999-BUL-36, which 
has t^ ~ 30 days. However, as we describe in Appendix the parameters we derive are reasonable: the 
asymmetry is quite small, and only detectable due to the excellent data quality of the event. Since both 
models fit the data equally well, we conclude that we cannot reliably distinguish between them, although 
we favor the parallax interpretation based on the fact that the planetary fit requires a special geometry and 
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a parallax signal must be present at some level in all light curves due to the motion of the earth around 
the sun. We therefore conclude that we cannot robustly detect a planet from an asymmetry that is equally 
well-fit by parallax. This in turn implies that all planetary perturbations consistent with such an overall 
asymmetry should be ignored in the efficiency calculation for all events. Although we have not done this, 
we have performed simulations which demonstrate that by not doing so, we overestimate our efficiencies by 
only a few percent, which is small compared to our statistical uncertainties. The parallax and planetary fits 
to OGLE-1999-BUL-36, as well as a detailed account of these simulations are presented in Appendix 

The light curve of MACHO 99-BLG-18 displays a ~ 15 day anomaly of amplitude ~ 2%. Such an 
anomaly is longer than that expected from planets with q < 0.01, and we therefore systematically explored 
binary-lens fits with 0.01 < <? < 1. This uncovered a fit with q ^ 0.2 that is favored over the best-fit 
planet {q = 0.01, d = 0.8) by Ax^ = 22. Clearly we cannot claim detection of a planet when a roughly 
equal-mass binary model provides a substantially better fit. However, since Ax^ = 22 is below our normal 
threshold (Ax^ = 60), we must estimate the probability that in excluding MACHO 99-BLG-18 from 
the analysis, we have inadvertently thrown out a real planetary detection. Naively, this probability is 
exp(— Ax^/2) ~ 10~^, but we have already seen that unknown systematic effects generate a whole range of 
planet-like perturbations at the Ax^ ^ 50 level. An upper limit to the probability that a planetary light 
curve has been corrupted to look like an equal mass binary can be estimated directly from the data. It 
is P < /apriori/22 where /22 ~ 20% is the fraction of events with Ax^ < —22, and /apriori is the a priori 
probability that the event contains a planet that is being corrupted by systematic effects into a. q > 0.01 
binary, rather than a true q > 0.01 binary. This last quantity is unknown, but since we detect of order 10 
other binaries and no other planets, /apriori is certainly less than 50%. Thus P < 10%. This probability is 
smaller than the statistical errors on our resultant limit on planetary companions from the entire sample 
of events. Thus, excluding MACHO 99-BLG-18 as a binary causes us to overestimate our sensitivity to 
planets, but by an amount that is smaller than our statistical errors. 

Thus, out of an original sample of 43 events, we are left with 42 events (rejecting MACHO 99-BLG-18), 
and no viable planet candidates. Given this lack of detections, we can use the individual event detection 
efficiencies to determine a statistical upper limit to the fraction of lenses with a companion in the range 
of d, q parameter space that we explore. 



6.4. Detection Efficiencies 

The detection efficiency ei{d,q) is the probability that a companion with mass ratio q and projected 
separation d would produce a detectable deviation (in the sense of Ax^ < — Axthrcsii) the observed light 
curve of event i. Figure || shows ei{d, q) for our fiducial threshold Axtj^cgi, = 60 and all our events in the 
parameter range we searched for companions, 0.1 < d < 10 and 10^"* < 9 < 10^^. 

We have plotted e as a function of logd, which clearly reveals the d ^ d^^ symmetry inherent in 
planetary perturbations ( |Griest fc Safizadeh 1998 ; Dominik 1999b ). For low magnification events {uq > 0.1), 



the efficiency exhibits a "two-pronged" structure as a function of d, such that the efficiency has two distinct 
maxima, one at max < 1 and one at d^maxj and a local minimum at d = 1. The approximate locations of 
these maxima can be found by determining the separations at which the perturbation due to the planetary 
caustic occurs at the peak of the light curve. 



.±1 1 1 



2"oT2\/"o + 4. (22) 
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For planetary separations dj^max < d < rf7maxi ^^"^ caustics produced by the companion are within a radius 
uq of the primary lens, and are thus not well probed by the event. For high magnification events, e is 
maximized near d — 1. This is not only a consequence of equation (p2|), but also because the central caustic 
is being probed by the event. As expected, the detection efficiency to companions with any q and d <^ 0.2 
or d ^ 5 is negligible for nearly all events. 

Of the 43 events, 13 have very little detection efficiency: for these events, e{d,q) is larger than 5% 
for only the most massive companions, and never gets larger that 25%. For the most part, these low 
efficiencies are due to poorly constrained uq. Eight events, notably all high-magnification events, have 
excellent sensitivity to companions and exhibit e(d, q) > 95% for a substantial region in the (c?, q) plane. 
Our resultant upper limits on small mass ratio q <^ 10""^ companions (§||) are dominated by these 8 events. 
For the remainder of the events, the efficiency is substantial 25%) for some regions of parameter space. 
These events contribute significantly to the upper limits for mass ratios q ^ 10"'^. 

In Figure ^ we show the efficiency averaged over the lensing zone (where the detection efficiency is the 
highest), 

^1.6 

eLzAl) = / ^^id,q) dd , (23) 

"'0.6 

as a function of the logarithm of the mass ratio. For a model in which companions have projected 
separations d distributed uniformly in the lensing zone, ei,z,i{q) is the probability that a planet of mass ratio 
q would have been detected in Hght curve i. Also shown is eLz.i for a detection threshold of Axthrcsh ~ ^O^^- 
For this more conservative threshold, the efficiencies are 5 — 40% lower, though the threshold level is most 
important where the efficiency is smallest. 



Finite Source EfTects 



The results in §6.3 and §3.4 were derived under the implicit assumption that the source stars of the 



microlensing events could be treated as point-like. Numerous authors have discussed the effect of the finite 



size of the source on the deviation from the PSPL curve caused by planetary companions ( Bennett fc 
Rhie 1996| ; paudi fc Gould 1997| priest fc Safizadeh 1998| ; |Gaudi fc Sackett 2000| ; [Vermaak 2000| ). Finite 



sources smooth out the discontinuous jumps in magnification that occur when the source crosses a caustic 
curve, and generally lower the amplitude but increase the duration of planetary perturbation. Finite 
sources also increase the area of influence of the planet in the Einstein ring. Thus finite sources have a 
competing influence on the detection efficiency: significant point source deviations can be suppressed below 
the detection threshold, while trajectories for which the limb of the source grazes a high-magnification area 
can give rise to detectable perturbations when none would have occurred for a point-source. Which effect 
dominates depends on many factors, including the size of the source relative to the regions of significant 
deviation from the single-lens form, the photometric precision, and the sampling rate. For large sources 



and small mass ratios, finite source effects can significantly alter the detection efficiency (Gaudi & Sackett 



200C). Since in principle the results presented in §§ |6.3| and 3.4 could be seriously compromised by ignoring 
these effects, we evaluate the magnitude of the finite source effect explicitly. 

In order to access the magnitude of the finite source effect, we must estimate the angular radius of the 
source in units of 9^, 



0.02 



6^asy V 12.5 kms-^kpc 



tE 



40 days 



(24) 
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where 0* = 6 /las for a clump giant at 8 kpc. For deviations arising from the planetary caustic, finite source 
effects become important when 9^ is of order or smaller the planetary Einstein ring radius, 6p, i.e, when 

P* <; (Planetary Caustics). (25) 



The size of the central caustic is Uc ^ qdid — f ) ^ ( Gricst fc Safizadeh 1998 ). Thus finite sources will affect 



the magnification due to the central caustic when > qd{d — 1)^^. However, in order for the central 
caustic to be probed at all, the event must have an impact parameter wq < Uc. Thus finite source will affect 
deviations arising from the central caustic if 

P* ^ ""0 (Central Caustic). (26) 

The difficulty in assessing the effect of finite sources on the detection efficiency lies not in evaluating 
the effect for a given p*, but rather in determining the appropriate for a given event. This is clear from 
equation (p^: of the three parameters that determine p», is known from the PSPL fit, 0, can be estimated 
based on the color and magnitude of the source, but ^roi is unknown. Gaudi fc Sackett (2000| ) suggested 



several possible methods of dealing with this difficulty. The simplest is to assume for all events a proper 
motion equal to the mean proper motion (/irci), adopting the given by equation (|2j) with /Xroi = (A*rci)- A 
more accurate, but also more complicated and time-consuming, method is to integrate over a distribution of 
/Ltrei given by a Galactic model; this would imply calculation of finite source effects for many different values 
of yO*. Here we adopt the first approach and determine assuming /irci — (Pici) — f2.5 kms^^kpc"^. This 
value of (/iroi) is a factor of two lower than the expected mean relative proper motion for all lenses toward 



the bulge (Han & Gould 1995), and refiects the fact that our median t-^ is a factor of two larger than the 
median of all microlensing events toward the bulge and our belief that the larger time scales reflect the 
fact that we are preferentially selecting slower (rather than more massive or closer) lenses. We justify this 
assertion in To the extent that the masses and distances of the lenses contribute somewhat to this larger 
median time scale, our adopted value of (/Zrei) is likely an overestimate. Therefore, the resulting values of 
are likely overestimates, so that we are conservatively computing upper limits to the effect of the finite 
source sizes on our conclusions. 



7.1. Estimating the Source Sizes 

The angular size 0, of a given source can be estimated from its dereddened {V — I)q color and 
magnitude Iq. From the PSPL fits to the / and V photometry, we know the / and V fiuxes of the sources 
in instrumental units (see §||). We assume that the dereddened color {V — /)ci,o and magnitude /ci,o of the 



clump is invariant for all our fields, adopting the determination by Paczyhski et al. (199£), 



(1^- /)ci,o = f.f 14 ±0.003, /ci,o = 14.43 ± 0.02. (27) 

We form instrumental color-magnitude diagrams (CMDs) for each of our fields, and determine the position 
of the clump in instrumental units by finding the local maximum in the density of sources. The difference 
between this position and the intrinsic position [eq. ([27|)] gives the offsets A{V — I) and A/ for all the stars 
in the field (except foreground stars, which have less reddening than the calibrating clump). Note that 
these offsets include both the calibration from instrumental to true fluxes, and also the correction for the 
mean reddening of the field. Thus we do not assume a constant redenning law from field to field. We apply 
these offsets to the instrumental {V — I) and / of our source stars, finally arriving at the {V — I)o and Iq for 
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all our sources. These are shown in Figure and listed in Table |. The error bars on {V — I)q and Iq are 
those derived from modeling uncertainties; we estimate there to be an additional calibration error of ~ 5% 
in both {V — I)o and Iq based on the typical dispersion of the clump. Note that the majority (~ 70%) of 
our sources are giants. 

Using these colors and magnitudes, the angular size of the sources are derived using a modified version 



of the empirical color-surface brightness relation derived by van Belle (1999 ) and given in Albrow et al. 
(2000a ). The resulting 9^, for all of our sources are shown in Table ^. The average uncertainty in 6^ 
is 0(20%), combining the uncertainty in the color and magnitude of the source due to both modeling 



and calibration uncertainty and the uncertainty in the underlying van Belle (1999) relation. We do not 



determine the uncertainty on 0* for individual sources because the uncertainty in p,, (the parameter in 
which we are primarily interested) is dominated by the uncertainty in (/irci)- Seven of our events have 
insufficient F-band data to determine the instrumental {V — I) of the source. For these events, we assumed 
the source to have the median (V — I) of all sources in the field with similar / magnitudes. (For these 
events, we do not quote uncertainties on {V — /)o-) Finally, four events had either no y-band data at all, 
or the position of the clump was impossible to determine from the CMD of the field. For these events, we 
simply adopt the conservative assumption that the sources are clump giants, with 6** = 6 /xas. 

These estimates of 9^ are used to determine under the assumption that all events have the same 
relative proper motion (/irci) = 12.5 kms^^kpc^^; these values of are listed in Table ^. For two events, 
MACHO 98-BLG-35 and OGLE-1999-BUL-35, the value of p, estimated in this way is larger than the 
fitted uo of the event. In both cases, the derived values of are ruled out by the fact that, despite dense 
coverage at the peak, no deviations from the PSPL form are seen, as would be expected if p* > uq and 
the lens was resolving the source ( iGould 1994 iNemiroff fc Wickramasinghe 1994| ; |Witt fc Mao 1994)) . For 
these two events, we therefore assume that = up. In Figure |ll|, we plot up versus for all our events, 
along with the boundaries where finite source effects become important for both the planetary and central 
caustics [eqs. (^ ) and (26)]. For the majority of our events, finite source effects should not alter the results 
for companions with q ^ 10^"^, whereas a large fraction of our events should be affected for q ~ 10^^. 



7.2. Incorporating Finite Sources 



In order to incorporate finite sources into the analysis, we repeat the algorithm presented in §S.l for 



all events, but fit the events to binary-lens light curves that include the effect of the finite size of the 
source. Evaluating the finite-source binary-lens magnification for the specific value of p^ determined for 
each event is not computationally feasible, as finite source magnifications are quite time consuming to 
calculate. We therefore adopt a procedure similar to that described in |6.2.4 : interpolation between a 
grid of finite-source binary- lens magnification maps. We choose the same grid spacing and size for {d,q), 
namely 10"** < q < 10^^ at equal intervals of 0.25 in logg, and 0.1 < c? < 10 at c? = 0.1, 0.2, .... 1.0 and their 
inverses. For each of these {d,q) pairs, we create finite-source magnification maps for 10^** < p* < 10~^ 
at intervals of 1/3 dex in logp*. These maps have same the extent and sampling in the source plane as 



the point source maps (see §8.2.4). We evaluate the finite-source magnification using the Stokes method 
of integrating over the boundary of the images ( [Kayser fc Schramm 1988t [Gould fc Gaucherel 199^ ). Our 
assumption of uniform sources overestimates the size of the finite source effect relative to limb-darkened 
sources, and thus is conservative. The grid value of p* closest to the value estimated for each source is used 
to calculate the detection efficiency for that event. We have repeated this process for the next-closest value 
of p, in the grid for all events, and find that there is no appreciable difference in the conclusions. 
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7.3. Effect of Source Size on Detection Efficiencies 

The distribution of ^Xmin ^"-"^ finite-source binary-lens fits is shown in Figure ^ along with the 
distribution for the point-source binary-lens fits. For the most part, the two distributions are quite similar. 
The significance of the best-fit binary-lens model has increased in some cases (e.g. MACHO 98-BLG-35), 
but all of the events that fall below our detection threshold {Axf^^i^^ > —60) in the point-source case also fall 
below this threshold in the finite source case. We recover the same two anomalies in MACHO 99-BLG-18 



and OGLE-1999-BUL-36, but no others. As argued in §6.3, these two anomalies have explanations other 
than planetary microlensing for their behavior. Thus our conclusions are unchanged: out of a sample of 43 
events, we find no viable planet candidates. 

The resulting finite-source lensing zone detection efficiencies [eq. (^3|)] are shown in Figure |^ along with 



the corresponding point-source efficiencies. We find, in agreement with the expectations in §7.1, that the 
difference between the point-source and finite-source efficiencies for mass ratios q ^ 10^"^ is negligible for 
nearly all events, with the exception of a few events with very large sources (p* ~ 0.1). Finite source effects 
begin to become appreciable for q <^ 10~^. For q = 10~*, the finite-source detection efficiency is markedly 
smaller than the point-source efficiency for large sources. The finite size of the sources has no appreciable 
effect on the detection efficiencies for those mass ratios where we have significant constraints {q > 10^^), 
and conversely for those mass ratios where finite source effects are appreciable we have no interesting 
constraints. Therefore we conclude that, for this sample of events, finite source effects are negligible. 



8. Upper Limits on Planetary Companions 

The fact that a large fraction of our final sample of 42 microlensing events has significant detection 
efficiencies to planetary companions — despite the fact that we have detected no viable planetary candidates 
in these events — suggests that the fraction of primary lenses with planetary companions in our range of 
sensitivity must be considerably smaller than unity. To quantify the exact limit implied by our data, we 
combine the individual event efficiencies {d, q) to obtain a statistical upper limit on the fraction of lenses 
with companions as a function of mass ratio q and projected separation d. 

Assume that a fraction f{d,q) of primary lenses have planets with parameters {d,q). Averaged 
over a large number of events, the probability that any single event would harbor such a planet is then 
also f{d,q). The probability that such a planet would be detected in event i is the detection efficiency, 
ei{d,q). Therefore, the probability that any given event has a planet that is detectable with these data is 
f{d,q)ei{d,q). The probability that a planet is not detected is 1 — f{d,q)ei{d,q). Thus the probability that 
a sample of N events would result in at least one detection is simply 

P{d, q) = l- nl, [1 - /(d, q)e,(d, q)] . (28) 

The 95% confidence level (c.l.) upper limit to f{d,q) implied by such a sample of events is found by setting 
P(d, q) — 0.05 and solving for /(d, q). Note that, in the limit of /e^ <C 1, equation ( ^ reduces to the naive 
formula, 

P{d, q)^l- exp[-iVexp(d, q)] iVcxp(d, q) = /{d, q) ^ e,;(d, q). (29) 

i 

We have, however, used the exact expression equation ( psf ) to compute excluded fractions f(d,q). 

In Figure ^ we show the 95% c.l. upper limit to f{d,q) as a function of d,q derived from our final 
sample of 42 events, assuming Axthj-csh — ^0 point sources. We conclude that < 28% of lenses have a 



- 26 - 



companion of mass ratio q ;> 10~^ and projected separation d ^ 1. The hypothesis that more than one- half 
of the primary lenses have a companion near d = 1 for the full range of mass ratios 10"^ < Q < 10^^ is 
excluded with 95% confidence. Also shown in Figure |l2| are cross sections of the (d, q) exclusion diagram 
(95% c.l. upper limits as a function of d) for three different mass ratios, namely q — 10^^, 10^^, and 
10^'*. For these cross sections, we also show the 95% c.l. upper limits derived assuming point sources and 
^Xthrosh — 100, and assuming finite sources and Axthrcsh — 60. Clearly finite source effects are negligible in 
regions where we have interesting constraints. 

In Figure ^ we show the 95% c.l. upper limit as a function of q for companions anywhere in the lensing 
zone 0.6 < d < 1.6, and anywhere in the "extended" lensing zone, 0.5 < d < 2.0. Statistically, less than 20% 
of primaries have a. q = 10^^ mass ratio companion in the lensing zone. For q = 10"'^ companions in the 
lensing zone, the upper limit is 45%. 



9. Converting to Planetary Mass and Orbital Separation 

The upper limits presented in §^ are the most direct, least model-dependent inferences we can draw 
from our data. Unfortunately, they are not the most illuminating, for several reasons. First, the nature of 
primaries around which we limit planets is not specified. Second, our results are quoted in terms of the two 
natural binary-lens parameters, the mass ratio of the system q and the instantaneous projected separation 
d of the companion, rather than the more common (and more interesting) parameterization of planetary 
mass rup and orbital separation a. 

Unfortunately, it is not possible to directly determine the mass of the primaries, and hence their nature, 
because the one observable parameter containing information about the lens, the event time scale ^e, is a 
degenerate combination of the mass, distance, and velocity of the lens [Eqs. |l| and Only model-dependent 
inferences about the nature of the primary lenses are possible. The majority of the microlensing events 
in our sample are likely to be due to bulge stars lensing other bulge stars (Kiraga fc Paczynski 1994|). 



Following Gould (2000| ), we adopt the bulge mass function as measured by Zoccali et al. (200C), and assume 



where r is the Galactocentric distance, and 
have Gaussian velocity distributions with dispersion a = 100 kms^^. This model gives typical parameters 
for bulge self- lensing events of (M) ~ 0.3 M©, (tTi-oi) = 40 ^as, and thus (^e) ^ 320 fias. For the relative 
proper motion, this model predicts (/Xrei) ~ 25 kms~^ kpc~^, and thus (^e) ~ 20 days, which is the median 



time scale found by OGLE for events toward the Galactic bulge (Udalski et al. 200C). Taken at face value, 
the fact that the median time scale of the events in our sample is a factor of two times larger implies that 
we are selecting a biased subset of lenses. From equations (^ and (||), this bias could be toward higher mass 
lenses, slower lenses (smaller /iroi) or closer lenses (larger n^ci), or any combination of these three factors. 
In fact, as demonstrated by [Gould (2000 ), the majority of the dispersion in the expected distribution of 



time scales arises from the dispersion in fi,-ci, not the dispersion in tTi-cI or M. This implies that we are, 
for the most part, preferentially selecting slower — rather than more massive or closer — lenses, justifying 
our assumption of (/iici) ^ 12.5 kms^^ kpc~^ for the estimates of in Thus the typical mass and 
lens-source relative parallax of the lenses in our sample is likely to be close to those of the complete sample 
of microlensing events. We therefore adopt (M) = 0.3 Mq and (TTrd) — 40 /xas, which for source stars at 
Ds ~ 8 kpc implies lens distances of -Dl ^ 6 kpc. In other words, the majority of our primary lenses are M 
dwarfs in the Galactic bulge. 



Some caveats must be noted. Kiraga fc Paczynski (1994) estimate that ~ 20% of events toward the 
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Galactic bulge are due to lensing of bulge stars by disk stars. Of the remaining ~ 80%, pould (2000| ) 
estimates that ~ 20% are due to remnants (white dwarfs, neutron stars, and black holes). Thus, we would 
expect ~ 60% of the events in our sample to be due to normal stars in the Galactic bulge. However, we have 
no idea which events comprise this 60%. Also, some fraction of the events in our sample are likely members 
of binary systems with separations that are either too small or (more often) too large to be distinguishable 
from single lenses. We have no way of determining which events these are, or even what fraction of our 
events are in such systems. Given our rather small sample of events and the uncertainties in the magnitude 
of these contaminations, we feel it is not appropriate at this stage to attempt to correct for these effects. 

The estimates of (Af) and (tTioi) adopted above imply (6'e) = 320 /xas and thus (rs) — 2 AU (for 
_Dl = 6 kpc) . We use these values to convert the upper limits derived in §H from dimensionlcss units to 
physical units, via the relations. 



Mjup, rp = (^^^ AU, (30) 
where rp is the analog of d (the instantaneous projected separation) in physical units. To convert from rp 



V0.003/ 



to the conventional three-dimensional separation a, we must convolve with the distribution function (Gould 
fc Loeb 199^ ), 

p(rp;a) = ^ 1-^ , (31) 
a \ I 

which is found by integrating over all random inclinations and orbital phases, assuming circular orbits. 
Thus the detection efficiency of each event i in the (a, mp) plane is, 

et{a,mp)^ / drpp{rp;a)e.i{rp,mp) (32) 
Jo 

These individual efficiencies ei{a,mp) can now be combined in the same manner as in §^ to derive 95% c.l. 
upper limits to the fraction /(a, rup) of events with companions as a function of the mass rup and separation 
a of the companion. 

In Figure |lj we show the 95% c.l. upper limit to f{a,mp) as a function of a and rUp, assuming 
^/^^thrcsh ^ ^^^'^ point sources. This figure is analogous to Figure ^2|, except that now our upper limits are 
in terms of the physical variables of the mass of the companion in Mjup and separation of the companion 
in AU, and we have identified our primaries as M-dwarfs in the Galactic bulge. In Figure ^ we show the 
95% c.l. upper limits to the fraction of lenses with planets in two ranges of orbital separations, (1.5 — 4) AU 
and (1 — 7) AU. Taking our inference about the nature of the primary lenses literally, we conclude that less 
than 33% of M-dwarfs in the Galactic bulge have Jupiter mass companions between 1.5 and 4 AU. Less 
than 45% have 3- Jupiter mass companions between 1 and 7 AU. These are the first significant limits on 
planetary companions to M-dwarfs, and are the primary result of this work. 



10. Discussion 

The majority of what we know about planetary companions has been gathered from radial velocity 
surveys of stars in the local neighborhood. However, these surveys have told us very little about planetary 
companions to M-dwarf primaries, as they have focused on F, G, and K-dwarf and have only recently begun 
surveying cooler stars. To date, the only M-dwarf with known planetary companions is Gliese 876 ( Marcy 
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et al. 1998; Marcy et al. 2001b). Our results therefore place interesting limits in an entirely new region of 
parameter space. However, this also means that the comparison between our results and those of radial 
velocity surveys is not entirely straightforward, as we are probing different primaries, and therefore different 
regimes of star, disk, and planet formation. Furthermore, our primaries are mostly old stars in the bulge, 



whereas those studied by radial velocity surveys are relatively young (Ford, Rasio, & Sills 1999; Gonzales 



1999; Santos, Israelian, & Mayor 2000). Finally, there is evidence that the host stars of local companions 
have super-solar metallicity (Gonzales 199£; Santos, Israelian, fc Mayor 2000| ), whereas stars in the Galactic 
bulge likely have solar to sub-solar metallicity. It is not at all clear how these differences between the parent 
samples we probe will affect the various proposed planet formation mechanisms. 

Rather than attempt to interpret our results in the context of these various parameters, which may or 
may not affect planetary formation, we simply make a direct comparison between our results and those of 
radial velocity surveys. In Figure |l6|; we show our 95% c.l. upper limits on the fraction of primaries with 
a companion as a function of the mass nip and orbital separation axis a of the companion, along with the 
measured mp sin i and a of those companions detected by radial velocity surveys. For the most part, radial 
velocity surveys are currently sensitive to companions of smaller a than is microlensing, although there is 
clearly some overlap. Also shown is the radial- velocity detection limit for a precision of 5 m s~^, a primary 
mass of 0.3 Mq (typical of our primaries), and a survey lifetime of 10 years. We also show the astrometric 
detection limit for 0.3 Mq primaries at 10 pc expected for SIM, which should achieve a precision of 10 /zas 
and have a survey lifetime of five years. 

The results from radial velocity surveys for companions indicate that / ~ 5% of local F, G, and 



K-dwarfs have companions between < a < 3 AU (Marcy, Cochran & Mayor 2000). It is interesting 
to ask how many more events we would need to monitor in order to limit the fraction of primaries with 
companions to 5% in the range of the separations to which we are sensitive. From equation (p9|), we find 
that, for small /, / cx N^J^. Given that our limits are / ~ 33%, we would require ~ 7 times more events of 
similar quality. This number could be significantly reduced if the quality of the alerts could be improved, 
i.e. if a larger fraction of events we monitor in the future were bright, high-magnification events. This will 



likely be possible with the next generation OGLE campaign (Udalski et al. 2000) 



11. Summary and Conclusion 

We have analyzed five years of PLANET photometry of microlensing events toward the Galactic bulge 
to search for planets. All of the 126 bulge microlensing events for which PLANET has acquired data over 
the last five years can be subdivided into three categories: events for which the data quality is too poor 
to determine the nature of the event, events that deviate from the single lens in a way not associated 
with planetary companions (roughly equal-mass binaries, parallax, finite source, binary source, etc.), and 
apparently normal point-source point-lens events (PSPL). We find no events in a possible fourth category: 
events that have short-duration deviations from the single lens light curve that are indicative of the presence 
of planetary companions to the primary microlenses. This indicates that Jupiter-mass companions to bulge 
stars with separations of a few AU are not typical. 

In order to justify and quantify this conclusion, we imposed strict event selection criterion, and derived 
a well-defined subset of 43 intensively monitored events which we carefully analyzed for the presence of 



companions. Using the method of Gaudi & Sackett (2000), we searched for the signatures of planetary 



companions in these events over a densely sampled, extensive region of parameter space. Specifically, we 
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searched for companions with mass ratios q from 10~^ — 10"^ and instantaneous projected separations d 
in units of the angular Einstein ring radius from 0.1 < d < 10. Based on an analysis of our photometric 
uncertainties for constant stars, we required that the difference in between the best-fit binary lens model 
and the best-fit single lens model be < —60 for a detection candidate. We found two such candidates, events 
MACHO 99-BLG-18 and OGLE-1999-BUL-36. Analysis of MACHO 99-BLG-18 revealed a significantly 
better fit with q ~ 0.2, and was eliminated from the sample. OGLE-1999-BUL-36 displays an overall 
asymmetry that is equally- well (in the sense of x^) explained by a low-amplitude parallax signal. Since 
we cannot reliably detect planets from global asymmetries, we explicitly discard this ambiguous anomaly. 
Thus we find no viable planetary candidates out of our original sample of 43 events. 

We then calculated the detection efficiency for our events in [d, q) space. Of our final sample of 
42 events (eliminating MACHO 99-BLG-18), 30 have substantial (> 25%) efiaciency for the detection of 
companions with q = 10~^ and separations in the lensing zone 0.6 < d < 1.6. Had all of the primary lenses 
harbored such companions, we should have detected a planet in at least ^ 7 of them. The fact that we 
detected no companions implies that this is not the case. By combining our efficiencies, we obtain statistical 
upper limits on the fraction of lenses with massive planets in the lensing zone. At the 95% confidence level, 
we find that < 25% of lenses can have a companion in the lensing zone with mass ratio q = 10~^ . 

Using a model of the mass function, spatial distribution, and velocity distribution of stars in the 
Galactic bulge, we infer that the majority of our lenses are likely due to M 0.3 Mq stars at 6 kpc, i.e. 
M dwarfs in the Galactic bulge. Using this assumption, we convert our upper limits from (g, d) space to 
mass-orbital separation space. We conclude that less than 33% of M-dwarfs in the Galactic bulge have 
Jupiter-mass companions between 1.5 and 4 AU, and less than 45% have 3- Jupiter mass companions 
between 1 and 7 AU. These are the first significant limits on planetary companions to M-dwarfs. 

We have also tested the robustness of our conclusions to various assumptions. The effect of the 

finite size of the source stars was estimated for each event using the color and magnitude of the source 
and assuming a mean relative proper motion of the lens. We find that the finite source efi'ect becomes 
important only for mass ratios q ^ 10~^, where our constraints on companions are already weak. We 
therefore conclude that finite source effects have a negligible effect on our results. We also tested the 
effect of changing our detection criterion from Ax^jji-cgh = 60 to Ax^j^j.^gj^ — 100. As expected, this lowers 
our sensitivity somewhat, and increases our upper limits by ^ 20%, but does not change our conclusions 
substantially. Finally, we have tested the effect of ignoring parallax asymmetries in the calculation of 
our detection eflaciencies, and find that this changes our limits by substantially less than our statistical 
uncertainties. 

We find that our median event time scale (^e = 40 days) is a factor of two larger than the median time 
scale for all events toward the Galactic bulge, a selection effect that arises from the manner in which we 
choose our targets. We argued that this primarily biases our events toward slower, rather than closer or 
more massive lenses. Therefore, our assertion of a typical lens mass of 0.3 Mq is justified. 

For the most part our upper limits are for planets with orbital separations that are larger than those 

currently probed by radial velocity techniques, since the orbital times are longer than the finite survey 
lifetimes. However, the smallest separations to which we are sensitive overlap with current radial velocity 
surveys, and as the radial velocity surveys continue, the degree of overlap will increase. Thus one will 
evcintually be able to compare the frequency of companions in the Galac;tic; bulge with that in the solar 
neighborhood. We estimate, however, that a sample ^ 7 times larger than that considered here would be 
needed to probe fractions as small as those being measured by radial velocity surveys (~ 5%), assuming 
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assuming that future microlcnsing observations are of similar quality to those analyzed here. If the number 
of alerts is increased substantially, however, more care could be taken to choose higher-quality (brighter, 
higher maximum magnification) events. This would considerably reduce the number of event needed to 
probe companion fractions of 5%. 

Our results have implications for theories of planet formation, as the orbital separations we probe may 
be closer to the sites of planet formation than the small separations at which radial-velocity companions 
are found, which may be reached via orbital migration. In any case, the limits described here provide 
fundamental constraints on the frequency and distribution of extrasolar planets orbiting the most common 
stars in our Galaxy. 
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Excluded Anomalous Events 



In §^ we rejected from the analysis 19 anomalous events which we asserted were not caused by 
planetary (i.e. small mass ratio binary) lenses. Here we list each of these events, and briefly justify why 
we believe their anomalies to be non-planetary in origin. For those events for which binary-lens fits are 
available in the literature, we will simply state the fitted mass ratio(s), and refer the reader to the paper; 
for a large fraction of these events, we rely on the analysis and binary-lens fits of [Alcock et al. (2000 ). One 
caveat should be noted. It is known (Dominik & Hirshfeld 1996; Dominik 1999a; Albrow et al. 1999b) that 
binary lens events, even extremely well sampled ones, often have degenerate solutions ( Afonso et al. 2000| ). 
This is due to intrinsic degeneracies in the binary lens equation (Dominik 1999b). Finding all of these 
degenerate solutions to an observed light curve is highly non-trivial, due to the extremely sharp variations 
in with respect to the canonical parameters, although several methods have been proposed to deal with 
this difficulty ( |Di Stcfano fc Mao 1996| ; |Di Stefano fc Perna 19"97| ; [Albrow et al. 1999b|) . It is therefore 



possible, as Alcock et al. (2000| ) allow, that not all solutions have been found and thus that some of the 
events they analyze actually have planetary solutions that they missed. Based simply on examination of the 
data we find this unlikely, since the deviations from the PSPL form are gross and long duration, contrary 
to what would be expected from a small mass ratio binary. 



For caustic-crossing binary-lens events for which the source is resolved, we can use the following 
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argument to place a lower limit on the mass ratio q. The maximum magnification obtained when a source 
of angular size 6** crosses a fold caustic is ( Schneider, Ehlers, fc Falco 1992| ) 



1/2 



whereas for a cusp caustic, 



(Al) 



(A2) 



Here Mr is a factor that describes the characteristic scale of the caustic. For caustics originating from binary 
lenses with small q, this scale is of order the planetary Einstein ring radius, Op. Due to possible blending, 
the observed maximum magnification, A^^^ „. , (or A"^„^ „k„), is a lower limit to the true magnification, and 
by combining eqs. (|Ai|), (^, and (24), we obtain an approximate lower limit on q for a fold caustic crossing: 
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and combining eqs. ([A^), (|l|), and (|2J), we obtain a similar relation for a cusp crossing: 
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where ts.obs is the observed (i.e. blended) time scale of the event, which is always a lower limit to the true 
time scale. Since gmin is proportional to tE,obs squared, while gmin is proportional to ^^^^x obs ^'^ ^^'^ fourth 
power and A'^^^ squared, the limits in Eqs. A4 and A3 hold even in the presence of blending. The 
smallest sources in the Galactic bulge have 6** ~ 1 /xas, and the dispersion in fi^oi for bulge-bulge lensing 
is a factor of ^ 2. Thus, an observed fold crossing with ^^axobs <^ almost certainly due to binary 
lens with mass ratio q > 0.01. A cusp crossing with A'^^^ ^j^^ ^ 40 is almost certainly due to a binary with 
q > 0.01. In general, for reasonably well-sampled events, a cusp approach can be easily distinguished by 
eye from caustic crossing events. For disk-disk lensing, for which /iici ~ 5 kms^^ kpc^^, somewhat smaller 
mass ratios are allowed; however such events are generally rare. 



MACHO 95-BLG-12 Both PLANET ( |Albrow et al. 199^ ), and MACHO/GMAN ( |Alcock et al. 2000D 
data show a smooth double-peaked event, with both peaks having comparable duration. This 

However, the 



morphology suggests a weak binary lens or binary source ( Griest fc Hu 1992 ). 
achromaticity of the event favors a binary-lens interpretation, and we find that a binary-source model 



provides a poor fit to the PLANET data. We cannot uniquely constrain a binary- lens fit, but Alcock 
et al. (2000| ) find a binary-lens fit with mass ratio q = 0.47. The fact that the peaks are of comparable 
duration precludes a small mass ratio binary-lens (i.e. planetary) model. 

MACHO 96-BLG-04 MACHO/GMAN data show two nearly equal-duration deviations separated by 



500 days (Alcock et al. 2000). Both deviations are separately well-fit by a standard PSPL model. 



suggesting a widely-separated binary-source or binary- lens ( Pi Stefano fc Mao 1996| ). Alcock et al 
(2000| ) find q = 0.88 for their binary-lens fit. Regardless of the interpretation, the PLANET data on 



this event would not have passed our second cut, due to insufficient data. 



MACHO 97-BLG-28 We find only one viable model that fits our data for this event (Albrow et al 



1999a), with q — 0.23. Alcock et al. (2000| ) find a similar binary-lens model fit for their dataset, with 
q = 0.21. 
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MACHO 97-BLG-41 Our data for this peculiar event is well fit by a rotating binary-lens model with 
mass ratio q = 0.34 ( Albrow et al. 2000a ). Bennett et al. (1999| ) favor the interpretation that this 
event is a planet orbiting a binary lens. Our data are clearly inconsistent with their particular fit, 
although this does not preclude the possibility that some fit of this nature would explain our data. 
Regardless of the interpretation, this event is rejected because of the presence of the binary. 

MACHO 98-BLG-6 This is a long-timescale (> 100 days) event which shows global deviations from the 
PSPL form indicative of parallax. 

MACHO 98-BLG-12 MACHO/GMAN data indicate that this event likely underwent four caustic 

crossings, with each pair of crossings separated by ^ 40 days ( Alcock et al. 200C ). The MACHO/GMAN 



data have poor coverage of the first set of caustic crossings, but constrain the amount of time that 
the source was between the second set of crossings to be ~ 3 days. Due to its short duration, one 
might suppose that the second set of crossings was due to a planetary caustic. However, the first 
set of caustic crossings, combined with the fact that the event exhibits a rise toward the second set 
of crossings, makes this interpretation impossible. Indeed, Alcock et al. (2000| ) find that the event is 
well fit by an intermediate-topology binary lens with q = 0.68. PLANET acquired a few data points 
immediately after the second crossing, and data immediately after the fourth crossing continuing 
until the end of the event. Due to the fact that the PLANET data did not probe any of the caustic 
structures, we find that our dataset is reasonably well fit by a PSPL model. However, our data alone 
fail our Suq/uq cut. 



MACHO 98-BLG-14 Both the MACHO/GMAN dataset ( [Alcock et al. 2000D and the PLANET dataset 
show a highly asymmetric light curve with a "shoulder" and then a peak. Such a morphology is 
indicative of a weakly-perturbed binary-lens event, and as such is prone to degeneracies. In fact 
Alcock et al. (2000 ) find two fits, one with mass ratio q = 0.09 and the other with q = 0.22. However, 
the event deviates from the PSPL form for a large fraction (~ 40%) of its apparent duration, making 
a planetary interpretation unlikely. We performed a systematic search of binary- lens fits to this event, 



using our data and the MACHO data. We recover the fits reported by (Alcock et al. 2000), along with 
a few other fits of similar significance. The best-fit binary with q < 0.01 is ruled out at the Ax^ = 50 
level. 

MACHO 98-BLG-16 MACHO/GMAN data show a short duration peak, followed by an abrupt rise and 
a plateau at magnification ^ 10 that lasts ^ 8 days. Following the plateau, the event returned to 
magnification ~ 2 ( Alcock et al. 200C ). Although the coverage is poorer, PLANET data qualitatively 
confirm this behavior. This morphology is consistent with a caustic-crossing binary lens event in 
which the short-time scale peak is due to a cusp approach, followed by a pair of fold caustic crossings 
with the usual intra-caustic plateau. MACHO / GMAN data near the peak of the first fold caustic 
crossing have A'^^^ ~ 20, and thus constrain the event to be non-planetary by equation (A3). 
Indeed, Alcock et al. (2000) find a binary-lens fit with q = 0.68. 



MACHO 98-BLG-42 Alcock et al. (2000) find q = 0.33. PLANET data cover the second half of the 
event, including the falling side of a second caustic crossing. Our data of the second crossing show no 
evidence of a cusp approach, favoring a pure fold caustic crossing. The data near the peak of this fold 
crossing have A'^ 



icf 
inax.obs 



40; therefore the event must be non- planetary in origin [eq.(A3)] 



OGLE-1998-BUL-28 This event displays a double-peaked structure indicative of a weak binary-lens or 
binary source. This is seen in both OGLE and PLANET data for the event. Using the combined 
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dataset, we find the best-fit binary-lens model has q = 0.34 and b = 0.42. Normalizing the errors to 
this model, the best model in the range q = 10^'' — 10~^ has Ax^ ^ 19, and thus is excluded. 

OGLE-1998-BUL-29 PLANET data for this high-magnification (Amax ~ 50) event show deviations from 
the PSPL form near the peak of the event that are indicative of source resolution effects. We find 
that a point-lens finite-source model fits the data quite well. In contrast, we find that the best-fit 
point-source planetary model in the range q = 10^^ — 10~^ is a significantly worse fit (Ax^ > 100). 

MACHO 99-BLG-8 Similar to MACHO 98-BLG-6, this long-time scale event shows severe parallax 
effects. We also find short time scale variability in the source. 

MACHO 99-BLG-22 Although the PLANET, MACHO, and OGLE data show no obvious anomalous 
behavior, our PSPL fit to the combined datasets yielded a time scale of t-E ^ 900 days, leading us 
to suspect parallax effects might be present. In fact, we find that a fit with parallax improves 
significantly, and results in a much more reasonable time scale. This interpretation is confirmed by 



the analysis of Mao et al. (2001). This event is excluded since our algorithm does not currently allow 



the search for planets atop other microlcnsing anomalies. 

MACHO 99-BLG-25 MACHO data for MACHO 99-BLG-25 show a clear deviation from PSPL at early 
times, in the form a smaller amplitude, but nearly equal duration peak occurring before PLANET 
began monitoring the event. The fact that both peaks are of similar duration suggests that this event 
is likely due to a binary source, and exclude the possibility that it is due to a planet. Our data only 
cover the rise and fall of the second peak and are perfectly consistent with a PSPL model. In fact, we 
find that this event does not have a significant planetary signal, nor does it have a large detection 
efficiency to planetary companions. Therefore excluding this event has no significant impact on our 
conclusions. 

MACHO 99-BLG-47 PLANET data show a departure from the PSPL form lasting ^ 3 days near the 
peak. Detailed analysis of this event shows that the deviation is caused by a small separation binary 



with q ^ 0.4 ( lAlbrow et al. 2001cD . 



MACHO 99-BLG-57 MACHO data show a large, long duration deviation from the PSPL form that 
is likely due to a binary-source or binary-lens. PLANET has very little data on this event, and so 
cannot confirm or clarify the nature of this anomaly. 

OGLE-1999-BUL-ll A caustic-crossing binary lens; PLANET data resolve the second crossing. The full 
dataset indicate a pair of pure fold caustic crossings. The second (fold) crossing has ^^^x obs ~ 



thus the event must be non-planetary in origin [cq.( A3)] 



OGLE-1999-BUL-23 We find only one viable fit to this event, with q = 0.39 ( |Albrow et al. 2001aD 



OGLE-1999-BUL-25 PLANET data show a large positive deviation lasting < 1 day superposed atop 
an extremely noisy light curve. As we see a sharp change in the slope of the light curve immediately 
after this deviation, we conclude that it is due a caustic crossing of some kind. It is not clear whether 
this deviation is due to a cusp or fold caustic crossing. We therefore conservatively assume that it is 
due to a cusp. The observed magnification at the peak of this deviation is A'^^^ ~ 40, and thus 
from equation (|A4|), the deviation cannot be due to planet. 



OGLE-1999-BUL-42 OGLE data indicate a double peaked structure to the light curve, which is likely 
due to an nearly equal-mass binary lens or a binary source. The PLANET data cover the rise and 
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fall of the second peak, and are consistent with a single lens model. Regardless of the nature of the 
anomaly, the light curve would not pass the cut on the uncertainty in uq. 



B. OGLE-1999-BUL-36 and Parallax Contamination 

Figure |l^ shows the PLANET and OGLE data for event OGLE-1999-BUL-36, binned into 1 day 
intervals, along with the best-fit single-lens, binary-lens, and parallax asymmetry models to the unbinned 
data. The difference in between the binary-lens and parallax models is 4 (with the binary-lens model 
giving the worse fit); both models are favored over the PSPL model by Ax^ ~ 80. The best-fit binary lens 
model has q = 0.0028, d = 0.60, and a = 1°. We also find fits for other mass ratios and separations that 



are nearly as good. The parallax asymmetry fit (see Gould, Miralda-Escude, & Bahcall f 994 for the exact 



form) yields a measurement of the asymmetry factor, n, given by. 



V, 



K — fl,^ — sinAsini/), (Bf) 

V 

where v — w(-Ds/^ls) is the transverse velocity of the lens projected on the observer plane, il^ — 27r yr^^, 
vq ~ 30 kms^^ is the speed of the Earth, and A is the angle between the source and Sun at the time of 
maximum magnification. In the case of OGLE-f 999-BUL-36, sin A - 0.6. We find k = 0.002f ± O.OOOf, 
which implies, 

= f43±7 kms"\ (B2) 

sin (p 

Combining this constraint with t-^,, we find an upper limit to the mass of the lens as a function of the 
distance to the lens, 

M<0.1Mq^—^, (B3) 

where x = Di^/Ds- Thus, if the parallax interpretation is correct, the lens must be closer to us than ^ 4 kpc 
in order to be above the hydrogen burning limit. 

The primary lesson learned from the analysis of OGLE-f 999-BUL-36 is that we cannot robustly detect 
planetary companions based on global asymmetries, since they cannot be distinguished from low-level 



parallax. However, when calculating our detection efficiencies (§6.4), we excluded all deviations that 
produced A^^ > 60, including asymmetries. Therefore, our efficiencies are overestimated. In order to 
estimate by how much, we choose a well-sampled, high-quality event, OGLE-I998-BUL-f 4, that contains 



data on both rising and falling sides. We repeat the algorithm in §6.1 to calculate the efficiency of this 
event but simultaneously fit for both the binary-lens magnification and parallax asymmetry. This procedure 
removes all detections based on asymmetry alone. In Figure ^ we show the detection efficiency both with 
and without excluding such detections. The difference is quite small, a few percent, because a very small 
range of angles produce deviations consistent with asymmetries. The majority of our events have sampling 
and photometric accuracy that is poorer than OGLE-f 998-BUL-f4, in which case they will be less sensitive 
to asymmetries. We therefore conclude that this effect is negligible. 
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Table 1: All Events from 1995-1999 with PLANET data. 
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Official Alert Name 


Abbreviated Name 


Status 


Notes 
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MACHO 
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Passed 
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Passed 
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Passed 




MACHO 
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MB97028 


Cut 
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MACHO 
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Passed 




MACHO 
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Passed 




MACHO 
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Cut 
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MACHO 
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Cut 
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MACHO 


97-BLG-59 
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Official Alert Name 


Abbreviated Name 


Status 


Notes 


EROS BLG-1998-2 


EB98002 


Passed 


- 


EROS BLG-1998-4 


EB98004 


Cut 


Insufficient Data 


MACHO 98-BLG-l 


MB98001 


Cut 


Insufficient Data 


MACHO 98-BLG-5 


MB98005 


Cut 


Insufficient Data 


MACHO 98-BLG-6 


MB98006 


Cut 


Parallax 


MACHO 98-BLG-12 


MB98005 


Cut 


Binary 


MACHO 98-BLG-13 


MB98013 


Passed 


- 


MACHO 98-BLG-14 


MB98013 


Cut 


Binary 


MACHO 98-BLG-16 


MB98016 


Cut 


Binary 


MACHO 98-BLG-l 7 


MB98017 


Cut 


Uncertainty in uq 


MACHO 98-BLG-18 


MB98018 


Cut 


Uncertainty in uo 


MACHO 98-BLG-19 


MB98019 


Cut 


Insufficient Data 


MACHO 98-BLG-22 


MB98022 


Cut 


Insufficient Data 


MACHO 98-BLG-24 


MB98024 


Cut 


Insufficient Data 


MACHO 98-BLG-26 


MB98026 


Passed 


- 


MACHO 98-BLG-27 


MB98027 


Cut 


Uncertainty in uo** 


MACHO 98-BLG-28 


MB98028 


Cut 


Insufficient Data 


MACHO 98-BLG-30 


MB98030 


Passed 


- 


MACHO 98-BLG-31 


MB98031 


Cut 


Insufficient Data 


MACHO 98-BLG-33 


MB98033 


Passed 


- 


MACHO 98-BLG-35 


MB98035 


Passed 


- 


MACHO 98-BLG-37 


MB98037 


Cut 


Uncertainty in uq 


MACHO 98-BLG-38 


MB98038 


Cut 


Uncertainty in uo 


MACHO 98-BLG-39 


MB98039 


Cut 


Uncertainty in uq 


MACHO 98-BLG-40 


MB98040 


Cut 


Uncertainty in uq 


MACHO 98-BLG-42 


MB98042 


Cut 


Binary^ 


*OGLE-1998-BUL-13 


OB98013 


Passed 


- 


*OGLE-1998-BUL-14 


OB98014 


Passed 


- 


*OGLE-1998-BUL-15 


OB98015 


Passed 


- 


*OGLE-1998-BUL-18 


OB98018 


Passed 




*OGLE-1998-BUL-20 


OB98020 


Cut 


Insufficient Data 


*OGLE-1998-BUL-21 


OB98021 


Passed 




*OGLE-1998-BUL-23 


OB98023 


Passed 




*OGLE-1998-BUL-25 


OB98025 


Passed 




*OGLE-1998-BUL-26 


OB98026 


Cut 


Uncertainty in uq 


*OGLE-1998-BUL-27 


OB98027 


Cut 


Uncertainty in uo 


*OGLE-1998-BUL-28 


OB98028 


Cut 


Binary 


*OGLE-1998-BUL-29 


OB98029 


Cut 


Finite Source 


*OGLE-1998-BUL-30 


OB98030 


Passed 
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Official Alert Name 


Abbreviated Name 


Status 


Notes 


EROS BLG-f999-f 


EB9900f 


Passed 


_ 


EROS BLG-f999-2 


EB99002 


Cut 


Insufficient Data 


*MACHO 99-BLG-6 


MB99006 


Passed 


_ 


*MACHO 99-BLG-8 


MB99008 


Cut 


Parallax 


*MACHO 99-BLG-ff 


MB99011 


Passed 


_ 


*MACHO 99-BLG-18 


MB99018 


Passed 


Deviation near peak*^ 


*MACHO 99-BLG-22 


MB99022 


Cut 


Parallax 


*MACHO 99-BLG-24 


MB99024 


Passed 


_ 


*MACHO 99-BLG-25 


MB99025 


Cut 


Binary Source? 


*MACHO 99-BLG-34 


MB99034 


Passed 


_ 


*MACHO 99-BLG-37 


MB99037 


Passed 





*MACHO 99-BLG-42 


MB99042 


Cut 


Insufficient Data 


*MACHO 99-BLG-45 


MB99045 


Cut 


Insufficient Data 


*MACHO 99-BLG-47 


MB99047 


Cut 


Binary 


*MACHO 99-BLG-57 


MB99057 


Cut 


Binary Lens/Binary Source? 


*OGLE-f998-BUL-5 


OB99005 


Passed 





*OGLE-f999-BUL-7 


OB99007 


Passed 





*OGLE-f999-BUL-8 


OB99008 


Passed 





*OGLE-1999-BUL-lf 


OB990ff 


Cut 


Binary 


*OGLE-f999-BUL-f3 


OB990f3 


Passed 




*OGLE-f999-BUL-f4 


OB990f4 


Cut 


Uncertainty in uq 


*OGLE-f999-BUL-f6 


OB990f6 


Passed 





*OGLE-f999-BUL-17 


OB990f7 


Cut 


Insufficient Data 


*OGLE-f999-BUL-f9 


OB990f9 


Cut 


Insufficient Data 


*OGLE-f999-BUL-22 


OB99022 


Passed 


_ 


*OGLE-f999-BUL-23 


OB99023 


Cut 


Binary^ 


*OGLE-f999-BUL-25 


OB99025 


Cut 


Binary 


*OGLE-1999-BUL-27 


OB99027 


Passed 




*OGLE-1999-BUL-33 


OB99033 


Passed 




*OGLE-1999-BUL-35 


OB99035 


Passed 




*OGLE-f999-BUL-36 


OB99036 


Passed 


Global Asymmetry'^ 


*OGLE-f999-BUL-39 


OB99039 


Passed 




*OGLE-f999-BUL-40 


OB99040 


Cut 


Insufficient Data 


*OGLE-1999-BUL-42 


OB99042 


Cut 


Uncertainty in uq^ 


*OGLE-f999-BUL-43 


OB99043 


Cut 


Insufficient Data 



Table I: Continued 



"Also shows evidence for binarity 

''The lightcurve of MACHO 99-BLG-18 has a small deviation near the peak of the event that is fit by a nearly equal mass binary 

lens. It is therefore excluded from the final event sample. 
'^The lightcurve of OGLE-1999-BUL-36 has a global asymmetry that is equally well-fit by a planetary model and a parallax 

asymmetry model. See §6.3. 
* MACHO and/or OGLE data included in the PSPL fit. 



References . — (l) |Albrow et al. (1998| ); (2) [Alcock et al. (2000D ; (3) |Albrow et al. (1999aD ; (4) |Albrow et al. (2000aD ; (5) |Albrow 
3t al. (2001a). 
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Table 2: Point-Source Point-Lens Fit Parameters for the final event sample. 
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Table 2: Continued 



*Indicates those events for which MACHO and/or OGLE data were included in the PSPL fit. 

''to=time of maximum magnification; t£;=Einstein ring radius crossing time; uo='naimm.ma impact parameter. 
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Table 3: Data Characteristics for the Final Event Sample. 
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Table 3: Continued 



^Indicates high-magnification (uo < 0.1; Amax > 10) events. 

"(Tnied- the median photometric error; Aimed- the median sampling interval. 
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Table 4: Source Characteristics for the Final Event Sample. 
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Table 4: Continued 



angular size of the source; p*=estimated angular size of the source in units of the angular Einstein ring radius of the lens. 
See §7.1. 

''InsufBcient V-band data to determine the color of the source; the source is assumed to have the typical {V — I) for its magnitude. 
''No CMD available, or CMD inconclusive. The source is assumed to be a clump giant. 
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Fig. 1. — (a) The percentage uncertainty Suq/uq in the impact parameter is plotted versus uq for all events that pass our data 
quantity cut. The uo for those events with 5uo/uo > 200% are plotted as triangles. The dotted line indicates our cut on the 
fra<;tional uncertainty, Suq/uq = 50%. Events with Suq/uq < 50% axe included in the final event sample and axe shown as solid 
circles, while events with Suq/uq > 50% are discarded (crosses and triangles), (b) The lines show the cumulative distribution 
of Mo for all events (solid) and those events that pass our cut (dashed). The histograms show the differential distributions of 
uo for all events (unshaded) and those events that pass our cut (shaded). The left axis refers to the cumulative distributions, 
while the right axis refers to the differential distributions, (c) The line shows the cumulative distribution of Suq/uq (bottom 
axis). The histogram shows the differential distribution (top axis). 
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Fig. 2. — PLANET data for the events that pass our selection criteria {§ |4|). The magnification, (F — Fs)/Fb, is plotted as a 
function of normalized time, r = (t — to)/tEi for the "cleaned" light curves, i.e. with seeing systematics removed and rescaled 
errors. See § H. 
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Fig. 3. — (a) The fitted impact parameter, uq, is plotted versus tlie fitted time scale, ts, for those events that pass all our cuts, 
(b) The line shows the cumulative distribution of tE (left axis), while the histogram shows the differential distribution (right 
axis), (c) The line shows the cumulative distribution of uq (bottom axis), while histogram shows the differential distribution 
(top axis) . The dotted line is for a uniform distribution in uo . 
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Fig. 4. (a) The histo gram shows the differential distribution of sampling intervals (in hours) for our final event sample (left 
axis). The solid line shows the cumulative distribution (right axis), (b) Same as (a), except in units of tE- The vertical dotted 
lines indicate the approximate minimum sampling rates necessary for detection of companions of the indicated mass ratios. 
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Fig. 5. — The median photometric error is plotted versus the median sampling interval in units of for our final event 
sample. The plus signs indicate high-magnification events («o < 0.1 or Amax > 10). The vertical dotted lines indicate the 

approximate minimum sampling rates necessary for detection of companions of the indicated mass ratios. The shaded box 
indicates approximately the median error and sampling needed to have significant sensitivity to 10~^ mass ratio companions. 
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Fig. 6. The vector positions in the source plane of the data points for event OGLE-1998-BUL-13, assuming a = 50°. We 
plot these for the best-fit ito as determined from the PSPL fit, = 0.30, as well as for the ±4(j bounds on uq. The arrow 
shows the direction of motion of the source with respect to the lens. The solid line connects the origin to the trajectory with 
uo = 0.30 at time t = to- Also shown axe contours of constant fractional deviation 5 from the PSPL magnification for a mass 
ratio q = 10~^ and projected separation of d = 1.11. The solid contours are <5 = oo, +5%, +1% (heaviest to lightest), while the 
dotted contours are 5 = —5%, —1% (heaviest to lightest). The solid black dots show the positions of the masses, the large dot 
is the primary lens, the small dot the secondary. 
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Fig. 7. — The shaded histogram shows the distribution of the difference in between the best-fit binary-lens model in 
the range q = 10~^ — 10"'*, and the point-source point-lens fit. Events with Ax^ < —40 are labelled. The dotted line is our 
adopted detection threshold, ^X^^^f.^-^ ~ ^''^ unshaded, bold histogram is the distribution of Ax^;jj found from a Monte 

Carlo analysis of constant light curves. See ^6.31 . (a) Binary-lens models in which the source is assumed to be point-like, (b) 
Binary lens models in which the source is assumed to have the dimensionless size p» given in Table kl 
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Fig. 8. — Black lines are contours of constant detection efficiency, e{d,q), shown for projected separations d between the 
primary and companion in units of the Einstein ring radius, of — 1 < log (d) < 1, and mass ratios between the primary and 
companion, q, of —2 > log(g) > —4. Contours mark e = l%(outer contour; dotted), 5%, 25%, 50%, 75%, and 95% (inner 
contour). Each panel is for a separate event; the abbreviated event name is indicated in each panel. The "wiggly" nature of 
the outer contours apparent in some events is an artifact of the (d, q) sampling and the plotting routine. Point sources have 
been assumed here. 
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Fig. 9. — Heavy solid lines axe point-source detection efficiencies averaged over the lensing zone (0.6 < d < 1.6) as a function 
of the mass ratio of the companion for a threshold of Ax^j^^^^j^ = 60. Dotted lines are for a threshold of Ax^^^J.^^^^ = 100. Dashed 
lines arc the lensing zone detection efficiencies for Ax^^^^^^^^ = 60 assuming a finite source of size in units of the angular 
Einstein ring radius. Each panel is for a separate event; the abbreviated event name and logp* are indicated. 
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Fig. 10. — The dereddened "clump-calibrated" color-magnitude diagram for the source stars. Crosses indicate sources for 
which separate lo and (V — I)o determination was possible; circles indicate events for which the {V — I)o was assumed to be 
that typical of stars with the same lo as measured for the event. The center of the clump is indicated with an arrow. Also 
shown is the CMD of a typical field (small dots). 
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Fig. 11. — The impact parameter «o versus the logarithm of the source size p* in units of the angular Einstein ring radius 6b . 
The dotted lines indicate the boundciries at which finite source effects become important for the detection of a companion of the 
indicated mass ratios via the planetary caustics; source sizes to the right of these boundaries significantly affect the amplitude 
and duration of the deviation caused by the planetary caustics. The dashed line indicates the boundary of the region at which 
finite source effects become important for the detection of a companion via the central caustic. 
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Fig. 12. — Left Panel: Exclusion contours (95% c.l.) for the fractions of primary lenses with a companion derived from our 
sample of 43 events, as a function of the mass ratio and projected separation of the companion. Solid black lines show exclusion 
contours for / =75%, 66%, 50%, 33% and 25% (outer to inner). The dotted (dashed) vertical lines indicate the boundaries of 
the lensing zone (extended lensing zone). Right Panels: Cross sections through the left panel, showing for three different mass 
ratios the upper limit to the fraction of lenses with a companion as a function of projected separation. The solid line is derived 
from the point-source efficiencies with a threshold of AXjj^j,gj,jj = 60. The dotted line is derived from the point-source efficiencies 
with a threshold of AyJ. . = 100. The dashed line is finite-source efficiencies with a threshold of AvJ. . = 60. The dotted 

^thresh ^^tnrosii 

vertical lines indicate the boundaries of the lensing zone 0.6 < d < 1.6. The dashed vertical lines indicate the extended lensing 
zone, 0.5 < d < 2. 
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Fig. 13. — Upper limits to the fraction of primary lenses with a companion as a function of the primary-companion mass 
ratio. Bold lines are for companions with projected separations anywhere in the lensing zone, 0.6 — 1.6 ve- Thinner lines are 
for projected separations in the extended lensing zone, 0.5 — 2 rg. 
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Fig. 14. — Same as Fig. [ij, except we have integrated over all possible orbital inclinations and phases to convert from 
projected separation to orbital separation, and assumed a primary mass of M = 0.3Mq and a primary Einstein ring radius of 
rE = 2 AU. 
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Fig. 15. — Upper limits to the fraction of primary lenses with a companion as a function of the companion mass. The bold 
lines are for companions with orbital separations 1.5 — 4 AU. The light lines are for orbital separations in the extended lensing 
zone, 1 — 7 AU. This figure is essentially identical to Fig. |l^, except we have integrated over all possible orbital inclinations and 
phases to convert from projected separation to orbital separation, and assumed a primary mass of M = O.SMq and a primary 
Einstein ring radius of te = 2 AU. 
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Fig. 16. — Our 95% c.l. upper limit to the fraction of M-dwarf primaries with a companion as a function of the mass rup 
and orbital separation a of the companion. The solid black lines show upper limit contours of 75%, 50%, and 25%. The points 
indicate the rup sin i and a of companions to stars (mostly G-dwarfs) in the local neighborhood detected by radial velocity 

surveys. Jupiter and Saturn arc marked with crosses. The dashed line shows the radial-velocity detection limit for a precision 
of 5 ms~^, a primary mass of 0.3 Mq, and a survey lifetime of 10 years. The dotted line is the astrometric detection limit for 
an accuracy of 10 /ias, a primary of mass 0.3 Mq at 10 pc, and a survey lifetime of 5 years. 
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Fig. 17. — Top Panel: Points show the magnification as a function of time for PLANET and OGLE data of event OGLE- 
1999-BUL-36, binned into 1 day intervals. The solid line shows the best-fit point-source point-lens (PSPL) model, the dotted 
line the best-fit parallcix asymmetry model, and the dashed line the best-fit binary model. Bottom Panel: The residuals from 
the best-fit PSPL model (in %) as a function of time. The dotted (dashed) line shows the deviation of the parallax asymmetry 
(binary-lens) model from the PSPL model. 
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Fig. 18. — Left Panel: Contours of constant detection efficiency, e, as a function of mass ratio and projected separation for 
event OGLE-1998-BUL-14. The contours are e = 5%, 25%, 50%, 75%, and 95% (outer to inner). The soUd contours are the 
efficiencies under the assumption that the parallax asymmetry parameter is zero, while the dotted contours are calculated with 
the asymmetry as a free parameter. Right Panels: Detection efficiencies as a function of projected separation for three mass 
ratios q. The solid lines are for no parallax asymmetry, and the dashed lines are with asymmetry as a free parameter. 
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